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1.  

e  eco:.. . 

^cs  of  the  syst.  &   too  fre 

system,  it  is  necessary  that  bhe  system  o]       with  s -- 

The  system  reliability  may  be  inc  tents 

which  form  the  structure  of  the  system.  This  -lied  upto  a 

certain  level  of  reliability,  t  . .  relial 

will  increase  the  cost  tremendoi .  rease  t 

reliability  for  a  system  composed  of  less  .  components  is  to  employ 

redundant  components.   In  case  of  parallel  redundan        r  of  c^ 

are  operating  simultaneously  out  of  which  only  one  cc.     at  is  r        for 

successful  operation  of  the  system.   In  case  of  -•--/■  er 

of  components  are  in  reserve  which  ..      switched  into  service  in  case  of 
component  failure. 

The  reliability  of  any  system  may  be  increased  by  employing  a 

^er  of  redundant  units.   Usually  there "e  restrictions  on  the  number  of 

units  th  .  .  ay  be  employed.  ?or  example,  increai  r  of  ei^. s 

results  into  increased  cost,  volume,  weight  etc.   The  resources  such  as  these 

being  limited  there  arises  a  rela  - 

two  types  of  problems  that  are  usually  encountered.   (l)  .  i  of 

system  reliability  subject  to  ir  and/or  nc ir  ■  i — .   (2)  Min- 

::ation  of  cost  of  the  system  while  achieving  a  c  Level  of 

reliability. 

The  purpose  of  this  report  is  to  obti 

by  the  discrete  ma:  Li    principle.   The  ;.         b  sin     r  multiple, 

.jar  or  nonlinear  cons  y  incre 

type  are  consider^.  I  • 


The  following  people  have  worked  on  these  probl i  with  linear  con- 
straints. Hees  and  Meerendonk  (6)  used  graphs  10  obtain  an  approximate 
solution.  The  lagrange  multiplier  technique  is  used  by     >witz  and 

McLean  (10)  and  Mine  (9).  Kettelle  (7)  used  dyn ic  progr ling  to  obt — . 

an  exact  solution  to  one  constraint  problem.   Belli  (3)  solved 

a  two  constraint  problem  by  dynamic  programming  using  lagrange  multiplier 
to  convert  the  problem  to  one  constraint  problem.   The  work  of  Kettelle 
was  extended  for  multiple  constraints  by  Proschan  and  Bray  .!_..   rhere  are 
number  of  defects  with  tneir  algorithm  such  as  very  few  cc;..  .-  mts  -— s  be 

handled,  due  to  approximation  made  for  the  computation  of  reli. „',  the 

optimum  solution  is  not  guaranteed  etc.   In  addition,  the  defect  of  the 
computational  efficiency  is  mentioned  in  their  paper. 

Even  though  the  solution  to  a  probl   a   y  appear  to  be 
practical  in  the  sense  -chat  tl     ;ic  quantities  do  not 
exceed  the  program  limitations,  it  nay  turn  out  to  be 
impractical  because  the  dominatin   squences  developed 
exceed  computer  capacity.   It  is  also  possible  that  th 
computer  time  required  y  be  excessive. 

Tillman  (13)  using  integer  pro^        solved  the  pj     is  of  maxi. 
systems  reliability  subject  to  mu   i        rable  constraints.   One  of  t 
problems  he  solved  is  similar  to  Example  5  of  this  report. 

Kettelle  (7)  vorked  with  a  nonlinear  constr problem  of 

a  linear  cost  function  while  achieving  a  i 'tain        level  -     ty. 

This  problem  is  the  same  as  Example  3  i  -  report.   - L3) 

applied  integer  progr     ,  for >n  of  a  . 

while  achieving  a  minimum  level  of  reliability  subject  to 

A  similar  problem  of  I     zing  system  re b ibject  to  noi ear  con- 
straints is  solved  i: 6. 


. 

problem  was  made 
lich  is  i.r.e  ex 

.   section  2  of      report  the  I 
review  of  systems  reliability  -..   .  ise 

in  system  r     Llity  optimis  basic  al- 

gorithm or  a  discrete  \  ie  i — ximui        Le  and  th« 

of  the  basic  algc.  L( i  ..re  give r -\...ce 

relation  for  a  one  dimensional  process  is  — •_" 

Section  k   deals  with  single  constraint  pre 

systems  reliability  subject  tc     ht  constr _s  giver.  Li         1.  Ti 

same  problem,  when  the  maximum  elements  at        -  is  rest] 

tc  some  value,  is  solvea  in  pie  2.  Both  1 -  L  with  r 

constraints.   In  Example  3,  a  ! ■  .on  is  . 

achieving  a  certain  minimum  level  of  reiia  with 

maximization  of  system  reliability  subject  tc  a  . linear  constraint. 

Optimization  of  problems  invc 

version  of  the  i     im  principle  is  c ^_  in  se      5«   

developed  in  this  section  is  \a be  s  Lve  syst< i  reli  —  lity  a\ 

problems  with  multiple  linear  and  nonlS     constr ts  in  section  6.   In 

j  5»  the  maximisation  o^"  I ility bt, 

cost  and  volume  coi  _e  6  deal- 

;  nonlinear  constraints  are  cons ' 

it  must  be  r.oiLu  I — t  the  number  ol  dur 5  is,  in 

reality,  positive  integer,  however,  -  luous 

variable.   Therefore  tl     Lutions  sj  iables, 

ley  are  n                                      rs.  It  is                         ations 
oot« 3  by  rou ;  i       still  *~_  — . —   ]  ..-  constr i, 


condition  of  convexity  or  concavity  is  a;    »d  for  insuri 

optimum.   If  these  tvo  as  .  searc. .  _   method  rna^  be 

used,  starting  with  the  solutj  in  th 

report.   Zero  one  integer  proga        armulatioi     n  in  i  lor.  7  may  be 

used  for  round  off. 


2.   . 

^ity  i~ 
successfully. 

.j vice  performing  its  purpose  u  I 
and  under  the  oper.  sred.  -.  of   reliabild 

is  ■  references  [l] 

In  practical  r<  -    Lity  - 

together  in  some  desil  -  reliability  . 

the  relationship  between  the  c.        of  as;  :t  on 

performance  of  -chat  system.   '_.  nship  c     -3  of  twc  first, 

the  output  of  the  individual  compone]        .  .    .   functional  - 

ship  among  ail  components  working  tof •  as  a  syste: 

result  of  the  operation  as  a  success  or  failure  is  cc. 

variation  in  the  quality  of  the  perfor ce  in  between  these  two  states. 

The  output  of  the  system  is  governed  by  the  fun     il  rela 
among  different  components.   The  sy~     Duld  be  desi  .  for  its 

operation  all  components  must  1  lly.         oyster.. 

_ed  a  series  system.      N-stage  series  system  is Pig.  1. 

other  hand  the  system  could  be  de:  ure  of 

any  one  of  the  components,  another  is  i  .  \,o   take  over 

rati  a;    refore,  the  syi  -'..:..".. 

for  successful  operation  at  least  one  :  form  si      rtorily. 

..  b.   system  configuration  is  calle.  -  _ 

system  is  shewn  in  Fig.  2.  'e  out 

N  simultaneously  or  ra  ssful 

operation  and  the  re:       N-l         s  serve  .-  re  sntly 
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d.  ■ 
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An   *■ 


described  "by  the  engineering  term         ic]  .   Most  of  the  prac — al 

systems  include  both  series  and  parallel  relat:  s; 

with  K  subsystems  in  series  -  si <syste:..  of  rec 

units  in  parallel  is  shown  in  Fig.  3. 

The  system  shown  in  Fig.  1  consists  of  a  number  of  components  comb: 

in  series  in  which  if  one  fails  the  system  fails.   For  ..   components  in 
series  the  system  reliability  is  given  . 

R  =  n  Rn  ,  (2-1)* 

s 

n=l 

where  R   is  the  reliability  of  the  nth  component  in  series  in  the  system. 

For  a  system  with  M  components  in  parallel  as  shown  in  Fig.  2,  — 
unreliability,  which  is  the  probability  that  all  M  components  will  fail, 

is  given  by 

M 

n  (1-r1)  . 


Therefore,  for  a  parallel  configuration  the  system  rel_ ity  is  given  by 

M 

R  =  1  -  n  (l-R1)  .  (2-2) 

If  the  reliability  of  each  component  of  the  parallel  system  is  equal  to  I 
then  equation  (2-2)  becomes 

-  l-(i--./-  .  (2-2a) 

3 

Now  if  the     N     stages   in  th  .   in   Fig.    3  h 

parallel  redundancies   at  the  nth 


Die  superscript     n     indicates   the   staj 

with  parentheses  or  brackets  such  as    (xn)     or    [13  l(^°     ;6n)j 
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K„  =  n  l-(l-Rnr  .  (2-3) 


NATURE  OF  THE  PROBLEM : 

Consider  the  example  of  the  three  s  syste.     .   —  - 

2  3 

ability  of  each  component  is  R"1"  =  D.8,  R  =  0.53     R  =  0.25.   rhis  gi'V 

the  total  system  reliability  as 

R  =  (0..     5)x(0.25) 
s 

=  0.1  . 

Now  if  we  add  one  redundant  unit  at  et s      e  system  reliability 

can  be  calculated  from  equation  (2-3),  which  is  given 

R  s  0.31?  . 

Thus  the  reliability  is  increased  by  a  factor  of  a  ...ere 

redundant  units  we  can  raise  the  level  of  the  r      lity  to  any  desired 

value.   Of  course  each  time  another  duplicate  component  1-  installed,  the 

total  weight,  cost  and  volume  increases.   Usually  t 

This  gives  rise  to  different  types  of  reliability  opti.     Lon  pre 

Type  1. 

The  problem  of  maximizing  the  rel: y  of  an  N-i       ries  s;  : 

with  redundancy  in  parallel  subject  to  multiple  lin     ad/or 

nonlinear  constraints  can  be  stated  as  follows : 

Maximize 

N 

r  =  n  ] 

subject  to 

t±{Bn)   <        i  -    2,  ....  r, 

where 


R     - 


.    - 


-      .e  reliabi 


)     =  the 


.    io"1)   =  relation  re  of  i  rce 

Cv..  .        

can  be  lir ■  — /or  qo r, 

r  =   ...  r  - , 

b .    =  the  total  amount 

i 

Often  there  is  a  restriction  on  the  iat  can  be 

used  at  any  one  stage.   In  such  a  case  the  aaditicnal  com 
written  as 

where  ia  is  the  maximum  numb  ,s   allowec  at  the 

Examples  1,  2,  k,   5  and  6  in  this  report  are  of  the  t        — ed 
above . 

Type  2. 

The  problem  of  option a  _i..^_.-  ot    active    fi 

a  certain  Ddnlmum  level  of  relii — ility  sd  as  J  s: 

Optica,  ze 

S  =   f(6ri), 

subject   to 

. 

n    Rn  >  R, 
n-l 

re 


f(9  )  =  the  objective  I  , 

R  =  the  level  of  2 

Example  3  is  similar  to  this  type.   It  should  be  

methods  used  in  this  report  can  be  applied  tc        3  of  this  type 

the  objective  function  is  non^ ear  an^-     -e  are  additional  lj -  - /or 

nonlinear  constraints  on  the  systi 


3.  :•■ 

-  process 

proces-  ...... 

is  transformed  at  each  stage   -      bo  i —  - 

-  ( 0-  ,0?,  . . . ,  6  ) ,  vh  . 

The  transformation  equation  at  nth  stage  may  be  writl , 

as  follows 

xn  =  r-,..--,  6n),         n-1,2,.,.,1,       -1) 

0 
x  =  a. 

The  optimization  problen  assc^_  .  -..~  a  seqi 

of  decision  variables  v\   n  =  1,  2,  .  .  .  ,  ..,  subject  to  const 

^(ej.ejj,  ...,  O  <  o,    b-1,2 :;,    ..-2) 

i  =  1,  2,  .  .  .  ,  r, 

which  makes  a  function  of  final  state  v   j   les 

s 

S  =  y  c.  :■:"'  c.  =  c     nt .  (5- 

,  *•-   x  i  i 

an  extremum  when  the  initial  co]      i  x  =  a  is 

procedure  for  solving  b m  optimiza  y  a  discr 

version  of  the  maximum  principle  is  to  ir.troc       s  di iional  a 

vector  z   and  a  Hamiltoni        wh 

relations 


c 

o 


- 


.. 


ev- 


er- 


cv 


T> 


••-  =   I  >  2, 


-1  ' 
3x. 

l 


_ 

J.  ,    w  i 

=  • 


t  o 


N 

-  c. 


4  > 


i  =  1,  2,  ...,s.    (3- 


If  the  opti:  .  vectoi  I 

function  3  an  optimum,  is 

(the  set  given  by  equation  (3—2)]  3  to  t 

(local)  extrer.:u=  with  respect  to  5   is 


56 


□ 


.". 


=  0,     a  >  1,  2,  3,  •  •  •  i  N. 


-.. 


If  6  is  at  a  boundary  of  the  set,  it  car.  be 

that  H   is  0-ocally)  an  ex 

This  is  the  basic  al  — 

above  formulation  of  the  problem  ve  s 
state  vard    a  are  fixed, 
final  conditions  of  th 

final  conditions  i ve 

ction  given  as 


16 


S  -       3.X?  . 

1  -  - 


Under  such  conditions  the  basic  algorithm  is  still  ap] ble,  ..*~z 

equation  (3-6)  is  changed  to 

z±   -  c±t  i  ■  1,  2,  .  .  .  ,  s, 

OPTIMIZATION  0?  ONE  DIMENSIONAL  PROCE    5  (k): 

If  a  multistage  decision  process  _d  for 

the  purpose  of  optimization  by  a  single  state  variable,  the  process  is 
called  a  one-dimensional  multistage  decision  process. 

For  one  dimensional  process,  there  is  only  one  state      ible  x.  , 
satisfying  the  performance  equation 

xj  =  T^xlJ"1;  6n),     n  =  1,  2,  .  .  .  ,  ...         [3-10) 

In  general,  the  objective  function  to  be  c      sd  is  sum  of  a  certain 

function  of  x.   and  6  over  all  ctajes  of  the  systei     -  as 

!  guv"1,  en)  . 

n=l 

The  optimization  problem  associated         .  a  pre. 

sequence  of  decision  variables   6'  ,  n  =  1,  2,  .  .  .  ,  ..   :o  -:  ig 


- 

L     G(x   i  e  ), 

n=l 


..  —  _ 

Xg  ■  x         +  G  ;   e  L,  2,    .    .    .    ,  N, 

■g-o. 

Therefore,   we   see   t  by 


B-  "         "      -.      .  [3-12) 

B-l 

Thus  the  prol 
sequence   of  6    ,   n  =  1,   2,    .    .    .    ,   N      is  to   be 

for  a  process   described  by  equati 

primary  state  variable  and     x       l  secon 

Tnen  the  Kami  It  on  i  an  function     .."'     defin  quatioi         -  ..  be 

written  as 


D    m    i    n-1         .-  ..-_  -  0-1 

..      =   z     T    (x.       ;    9    j   +  zg    [Xg       +  G   (x.       ;    8   ))    .  -13j 

According  to  equation  (2-5),  th  

vector  elements  z.   ana  z„  are  found  zc 


si 

n— i 

zi  =  —.t:  = on —  -  + — -— —  *2«     --■ --■ 


I*  ~  X  }   t  j   •   *   ■ 

Since  the  objective  function  is 


S=  I     Ci  Xi  =  X2  • 
1=1 


cx  =  0,  c2  =  1 


Thus  we  obtain 


z-L  =  0,  -15a) 


z2  =  1-  (3-15b) 


Combining  equation  (3-15b)  and  (3-15)  and  substitut  [3-1*0, 

gives 


zS  -  1,     n  -  1,  2,  ....  , 


and 


- 


1 


^ 


■ T « 


Combining  equations  (3-l6)  and  (3-13)i 


n  r,  /  n— 1       0  i  n-1 
--   -      (xx  ;  6  )  +  G       .  j ■  +  x   , 


r.  -  1 ,  2 ,  .  .  .  ,-•• 

According  to  equation  (3-7) i  i.e.,  th 

6  nay  be  found  where 


=  Zl       i?        +       "7 


Solving  this  equation  for  z. ,  we 


;  8  ) 


n 

Zi  " 

n 

-• 

- 


Substitution  of      Lon  [3-3 
relation 


<M 


8G  (x,   ;  6  )     3G  (x.  ;  8 


aen  aen+1  !;  r""' 


3T  (x"-x;  6")     3T  (xt,;  C"'A)       3xt" 

i 


n  ■  1,  2,  .  .  .  ,  N-1. 

Combining  equations  (3-15a)  and  (3-l8)  ^ives 


3G  (x?"1;  Bn) 

±-r= =  0.  (3-20) 

c/6 

Making  use  of  the  recurrence  relation,      ion  (2-19),  l  the 

performance  equations,  equation  (3-10)  and  equation  (3-20^  5  _  ..      of 

optimization  problems  associated  with  one-dimensional  processes  can  be 

solved.  For  processes  with  a  fixed  end  point  x. ,  con       z     =  0 

(equation  (3-15a))  or  equivalently ,  equation  (3-20)  is  deleted. 


-. 


EXAJ 


Lev 
sub-system  will  have   e  lei.      — 

problem  is   to  determine  the   o::  -  ^.-   i :h   ^_..- 

system  so  as  tc  Lze   the  total  -.  relia 

straint  that  the  total  m  in  the  rei — -  >s   ~..^ dl  gceed 

a  certain  fixed  amount. 

To  solve  this  problem  by  the   ^-~_.v~.-   .  i  us   c 

each   sub-system  as   a  s;.  _  let 


6     =  Number  of 


-• 


number  of  redun — r.cies  ^.  -l)i 

w     »  Wei    ...    of  one  redu.. 

x-  =  Amount  of  weight  left  ove     -;-  allocat     Lrst  ..  — 

xV  =  Lo ~  . 

n      -j . 
a  N  stage  system  can  be  represented 

n    n-1        . N          to 
x.  =  ..    -  w  (6  -1),      D  ■  1,  2 H,  v--- 

x 

x^  =  W, 
«?>0, 


x2  =  X2        [1-Cl-B  J  J,      n  =  1,  2,  .  .  .  , 

*2  -  °- 
)  is  constrained  to  positive  integers  only,  i.e., 

6n  >  1. 


The  objective  function  to  be  maxi:.'iized  is  c_  follows 


N  „n 

S  =  I     in      (l-(l-Rn)6  J  =  yC    .  (U-3) 

n=I 

N 
Xp  is  logarithm  of  total  system  reliability.  The  system  reliability  will 

be  maximized  if  we  maximize  the  logarithm  of  the  -•-" lity. 

Comparing  equations  (i+-l)  and  i.k-2)   with  the  equations  of  one  

process,  equations  (3-10)  and  (3-11),  ve  see  ti- 
me n~l  *n\         n~l    !i/,n  .  x  ,  v 
T(x   ;  6  )  =  x-   -  w  (8  -1),                       (4-J 

■L  a. 

G(xJ-1;  6n)  =  to   (l-(l-Rn)Qn)..  (U-5) 

Differentiating  equations  (4-1*)  and  (4-5)  with  respect  to  x? 


gives 


3T(x"-1;  6n) 

-=-wn, 


3Gn 


3T(x^1;  6n) 

=  1, 


a  n~1 
3X1 


- 


_- 


n-1 
3u(  x.       ;    I 

: : =     C     . 

..-1 


Let   us   define 


T,n       ,    „n 
u      =   1-R      . 


Substituting  equations   [k- 

one  dimensional  process,   equ 


) 


-    W  -    V 


,Tln-KUen+1 

(IT       ) 


i-(un+-)°  1- 


m+1 

-   1   = 


(un+.r 


or 


,  n   „   „n+l      u8  , 

ln\l-r~~    — 

i     n+1   „ 

en+I  = _ _ 

ji+1 
-  in  IT 

This  is  the  recurrence  relation  of  the  opti      iision.   It  must 
noted  that,  while  8  ,  n  =  1,  2,  .  .  .  ,  i:  are  in  reality  positive  -     rs, 
we  have  assumed  that  0  are  continuous  variables  in  obtaining  the  at 

rence  relation.  This  recurrence  relation  can  be  used  to  obtain  a  seqi of 

6  which  will  not  violate  the  weight  constraint.  The  trial  and  error  method, 
given  below  is  used  to  obtain  sequence  of  decision  variable. 

Step  1.   Assume  6  =  1  . 

2   3  N 

Step  2.   Compute  0  ,  0  ,  .  .  .  ,  0  from  equation  (4-10). 

N 
Step  3.   Compute  x  =  W  -  \     wn(0n-l)  . 

n=l 


N 
Cl 


x.,  will  be  in  either  one  of  the  fc       situations: 


than  zero,  (b)  equal  to  zero,  (c)  less        ro.   If  it  is  (a]  - 

go  to  step  kt   if  it  is  (b)  then  we  r< id  the  optimal  stage, 

to  step  6,  if  it  is  (c)  then  go  to  step  5. 

Step  k.      Increment  0  by  1  and  go  to  step  2. 

1 
Step  5-   Decrease  0  by  0.1  and  go  to  step  2  until  :. 

N 
zero,  when  x  is  greater  tb  i   ro  then  go  to  step  o. 

Step  6.  Solution  has  reached  the  optimal  s1  .  to  an 

integer  solution. 

A  computer  program  and  the  flow  <     m  for 
Appendix  A. 


A  five  stage  and 

program  .. 

number  of  stages. 

The  constants  as^-  le  ^. 

The  optimum  re.  led  is  i i  foU  . 

61  =  2.6000, 

62  =  3.8691, 

e3  =  5.2200, 

e*  =  3. 5606, 

65  =  3.053^. 
Since  6,n=l,2,  ...,5  should  he  positive 

*2  =  K 
e3=5, 

a5  =  3. 

should  he  noted  tta  I  ove    .'.  ;ures  icta 

stage;   to  obtain  the  number  of  re  -j  ve   si 

eac.'i    figure,    i.e.,    the   r.  . 
number  of  Ln 

The   total  wed  . 

which   is   equal  to  the   COB  □ is  gi 


2o 


Table  1.      Constants  assigned  for  5  sta^e  problem  . 

n R° £J V_ 

1  .90  8 

2  .75  9 

3  .65  6  101+ 
1+                               .60                              7 

5  .85  8 


The  cons^or.^3  assigned  for  a  ten 

The  op-  is  at 

.  DO, 
e2  =  1.9227, 
e3  =  1.691U, 
Qk  =  2.75J 
e5  =  3.UUU9, 
e°  =  1.7822, 

9T  =  2.9303, 

e8  =  1.660U, 

69  =  2. 3117, 
e10  =  1.91*96. 

Since,  9*,  n  =  1,  2,  .  .  .  ,  1  positj  -  -~ 

3 
9  =  2, 

•"-3. 

e5  =  i, 

a6  -  2, 
a7  =3. 


Table  2.   Constants  assigned  for  10  sta^e  problem. 

n R w 

1  .78  17 

2  .86  19 

3  .88  25 
k                             .80  8 

5  .70  8  200 

6  .90  16 

7  .65  21 

8  .9k  11 

9  .79  18 
10  .86 


e10  =  2. 

These  are  the  number  oj     ;s  at  i 

elements  at  each  stage  w 

The  total  weight  in  .....  9  is  1 

is  a  slack  of  2  units  ;ives 

of  0.7653. 


3C 


EXAMPLE  2:   MAXIMIZATION  07  THE  S 

CONSTRAINTS  WHERE  MAXIMUM  S  DF  REDUNDANCIES  

(LINEAR  CONSTRAINT). 

The  system  considered  is  the  same  as  in pie  1  out  —  

number  of  redundant  elements  at  any  stage         -  to  some  number  "..  , 

i.e.,  the  maximum  number  of  elements  at  any  stage  can  not  e;:_ ...  ••-  1. 

The  weight  constraint  remains  the  same  as  in  the  previous  problem.  The  con- 
straint on  6  can  be  written  as 

n 

G   <_  m  +  1  . 

The  performance  equations  for  this  problem  remain  the  same  as  in  — 
last  example;  thus,  we  obtain  the  same  recurrence  relation  of  the  optimum 

decision  variable.  The  only  change  is  the  procedure  of  sol-.      lis  eqi tion 

to  obtain  the  desired  result  in  the  light  of  the  new  co?.~\,~^-..\.   on  the dsion 

variable.  We  note  that  the  order  of  numbering  one  stages  is  arbitrary  and 

does  not  affect  the  optimal  decision.  We  wi__  le   use  of  this  . 

the  problem. 

Step  1.   Solve  the  problem  as  in  example  1,  without  any  considerat 

additional  constraint. 
Step  2.   Check  whether  6,n=l,2,  ...,.      less  t +  1. 

is  satisfied  then  the  new  restriction  does  not  a      cur  c 

decision  and  the  problem  is  solved.  But  if  for  some  ..      8 

is  greater  than  m  +  1,  then  go  to  step  3. 
Step  3.   If  there  are  p  stages  x'^r  which  6  is  less  than         ..  r v.  .  s 

the  stage  number  so  that  for  which  9  is  -^~^  tl     •-_" 

are  numbered  as  n  =  1,  2,  .  .  .  ,  p  and  tl  

0  is  greater  than  or  equal  to  "m+1" 

n  =  p+1,  p+2,  .  .  .  ,  :.  eu  /or  n  =  p-KL,  p+2,  .... 


Step  k.     Con 

H ,  total 

W  =  ..  - 

+1 

Step  5»  Go  to  step  1  with 
The  computer  program 
Appendix  I  . 

..RICAL      ..ES 

A  three  stage  and  a  ten  stage  prc~_< — i  are  -el/.'*..,  ...   restrj 
on  6  of  maximum  equal  to  five.      i_  ,  t) 
units  at  any  stage  is  restricted  to  four.   Hovev. 
developed  and  the  computer  program  are  for  zy. 

of  stages  and  the  number  of  reu cies  : 

It  should  be  noted  that  this  computer  j  i  .,  solv^ 

of  Example  1  by  putting  in  a  very  large  number  of  red  —  - 

constraint. 

...e  constants  assigned  for  a  —  -. 

The  optimum  redundancy  obtained  by  nu ric —  so]  a  be 

(see  Appendix  3,  Table  B3) 

61  =  2.37. 

62  =  5-0000, 

e3  =  5.0000. 


Table  3.   Constants  assigned  for  3  stage  problem. 


R  w  W         6 


1 

.95 

8 

2 

.80 

9 

3 

.65 

7 

75 


32 


,    n  =   1 ,    i  ,  ... 


-2, 


e2=>, 
a3  =5. 

It  is  i:.  :  tha  .   on    - 

nuaber  of  redundancy,   ve 

61  =  2.U00O, 

8     ■  U.0085, 

63  =  5-97      • 
As    6    ,   n  =  1,   2,    3  should  be  posi  .  .-s  we  c 

a1  «2, 

6     =   U, 

u3  =   f  . 

Since  the  above   figures  are  - ,   we   I 

subtract  one   froa  each  c.  e  to  c 

stage. 

note  that  when  the  Dumber 

of  four,  ve  obtain       J  wei 

system  reliability  of 

number  of  r<.     icy  total  weight  is 

reliability  of  0.99^1<   We  see  Li- 

is  reduj 


The  constants  assigned  for  the  ten  stage  sy;  ile  t. 

The  optimum  policy  obtained  by  numerical  solution  is  given  below  (see 
Table  B5).   It  should  be  noted  that  the  stage  numbers 

the  same  as  given  below.   This  is  because  coi      has  rearran^__  the  stage 
numbers  as  explained  in  step  3  above. 

e1  =  4.iooo, 
e2  =  3.2336, 
63  =  2.90^9, 

h 
e  =  5.0000, 

e5  =  3.9727, 

e  =  2.9041, 

eT  =  5.0000, 

e  =  2.5310, 

e9  =  3.9606, 

e10  =  3.2611. 

As  0  ,  n  =  1,  2,  .  .  .  ,  10  should  be  positive  integers,  we 


o1^, 


e2  =  3, 


a3  =  3, 


e5  =  4, 


Table   k.      Constants  assi.  or  10   . 


•• 


1 

• 

17 

2 

. 

19 

3 

.88 

25 

k 

.70 

li* 

5 

. 

15 

6 

.90 

16 

7 

.65 

8 

. 

9 

.79 

18 

D 

. 

18 

e6  =  3, 

e8  =  3, 

e10  =  3. 

As  indicated  "before,  to  obtain  the  number  of  redundancy  at  each  stage  ve  sub- 
tract one  from  each  of  the  above  figures. 

The  total  weight  in  redundant  units  is  ^66;  thus,  -/.ere  is  a  .      of 
7,  and  this  gives  the  optimum  reliability  of  0.97&2. 

Appendix  B,  Table  Bo  gives  the  result  of  the  same  problem  ' 6"  is 

unconstrained.  From  the  table  we  note  that  optimum  policy  reroains  ... 

as  before;  that  is,  the  additional  constraint  does  not  affect  the  optimum 

policy  for  this  case. 


a  a 

■ 

parallel. 

The  systt. 

its  efficJ  ..  .   

the  cost  of  the  syst< ...e  prc- 

elcc.u;::-j  while  ir.ai:. 

Consider  each  subsystem  as  a 

R  =  Mini     -vel  of  sved, 

r 

6*  -1=  Number  of  re - 

c'  =  Cost  per  redund~:;t  «_..  ;e, 

._■  logarit. 


x2  =  Cost  in  n  -s  up  to  and  inc] 

This  U     stage  system  car.  oe  represent<_ 


n    o-l 


_ 


xr»<*+ln{l-     T   }  ,     -  1,  2,  .  .  .  ,  .., 


- 


x?  =  0, 
1 


>  in  ( R ) , 


x  ,  =    x2    +  C  (  0  - 


~>  2, 


(fc-J 


x°  =  0. 


50 


3  should  be  a  positive  integer,  i.e., 


)n  >  1. 


We  see  that  the  objective  function  to  be  minimi is  given  by 

TT 

s  ■  I    4=    l    cn(en-D  =  jj  .  [k-ik) 

n=l      n=l 

Comparing  the  performance  equations  (4-12)  and  (4-13)  with  equation:;  of  - 
dimensional  processes  given  by  equations  (3-10)  and  (3-11),  we  fin 

T(x^-\en)  =  x^"1  +  to  {l  -  (i-Rn)en}  ,  (U-15) 


„/  n-l  nN    n,  n  ,  *  ,, 

G(x1   ;6  )  =  c  (6  -1)  4-16 

Taking  the  partial  derivatives  of  equations  (U-15)  and  (4-l6)  with  respi 

to  x         and  8     yields 


M*l      "'*>        -(l-**)6"   £n(l-R») 

l-d-H*)6" 


(4-.. 


— Vn — sl« 

8xl 


3G(x*-1;en) 
n 

38n 


3G(x"-1;en) 

~7^ — =0' 
3xi 


IT  =  :-.  ••  . 

one  dimensional  processes  given  b; 


-  -  "  ,un 


1-(U  )  1-     ) 


or 


>l,9n+1    ,   Irn+1  n+1 


=  {~ .   I" 


..  .n+1 
l-d/1*1)6  °      1- 


Comparing  this  with  the  rec  . 

the  same,  with  w  replaced  by   c  j      ,  ve  ol 


relation  as 


in  (1  ♦  [• 1)} 


This  is  a  recurrence  r 

noted  that,  while  9  ,  r.  =  -,  .1,  .  .  .  , 

have  ax,      that  o"  are  co:.  -  sur- 

rence  relation. 


This  recurrence  relation  can  be  uj     d  obtain  the  sequence 
8  which  vill  minimize  the  cost  in  redundant  elements  for  an 
level  of  reliability.  The  following  trial  and  error  method  is  used  to  obt 
the  sequence. 

Step  1.  Assume  9  -  1. 

2   3  N 

Step  2.   Compute  0  ,  6  ,  .  .  .  ,  6  from  equation  (4-21). 

Step  3.   Compute  x1  =  £  Jin  (l  -  (l-Rn)   )  . 

n=l 

x  will  be  either  (a)  greater  than  Zn(R)>  (b)  equal  to  Jln(. 

or  (c)  less  than  Jln(R).   If  it  is  (a)  go  to  step  5,  if  it  is 

(b)  go  to  step  T>  and  if  it  is  (c)  go  to  step  k. 

Step  k.      Increase  6  by  one  and  go  to  step  2. 

1  N 

Step  5.  Decrease  0  by  0.1  and  go  to  step  2  until  again  x*  is  less  than 

£n(R).   When  x.  is  less  £n(R)  go  to  step  6. 
Step  6.   Increase  0  by  0.1  and  compute  6,6,.  ...  ,6  from  equation 

(4-21). 
Step  7*  Approximate  0  ,  n  =  1,  2,  .  .  .  ,  N,  to  the  nearer  integer. 

A  computer  program  and  flow  diagram  for  the  IBM  lo20  is  given  in 
Appendix  C. 

NUMERICAL  EXAMPLES. 

A  five  stage  and  a  ten  stage  problem  is  solved.  However,  the  computer 
program  and  the  numerical  method  developed  i'j  for 
stages. 

The  constants  assigned  for  a  fiv 
The  optimum  redundancies  obtained  from  the  above  recurrence 
follows  (see  Appendix  C,  Table  C2): 


e3  =  3. 

e5  =  U.1501. 

As  6  ,  n  =  1,  .  .  .  ,5  she . 
62=5. 

e3  =  3, 

eU=3, 

These  figures  indicate  t..  To  c 

.oer  of  re at  units  at  e 

The  reliability  of  the  system  is  0*9023 
This  resu.-^  _:.  a 
...e  constanta 

ja  redundancies 
Table  C3). 

e1  =  2. 

62  =  -. 

65  - 


1*2 


Table  5.      Constants  assigned  for  5  sta~e  problem. 

Dn  n 
n « c 

1  .65  2.00 

2  .55  3.00 

3  .70  6.00                              .    DOG 

4  .75  7.00 

5  .60  i+.OO 


^e   u.      Constants   c  >r  10  i 

1 

2  .60 

3  .95 
.80 

5  .  . 

6  .90 

7  .80 
6 

9 
10  .65  2.- 


6°  =  1. 8U53, 

e7  =  2.5969, 

6°  =  5-2247, 

69  =  3.23^7, 

010  =  U.6136. 

As  6  ,  n  =  1,  2,  .  .  .  ,  10  has  to  "be  positive  integers,  ve  obtain 


e1- 

3, 

e2  = 

*. 

e3  = 

2, 

6^  = 

2, 

e5  = 

3, 

e6  = 

2, 

eT  = 

3, 

e6  = 

5, 

e9  = 

3, 

e10  = 

5. 

Tnese  figures  being  the  number  of  units  at  each  stage,  we     the  number  03 
redundant  units  by  subtracting  one  from  each  of  them. 

This  policy  gives  the  reliability  01*  the  Bystem  of  0.85^1  -     -  — 
requirement  of  O.65OO  and  the  total  minimum  co..  of  ..CO.  00. 




Lee,  etch  stage 

is  sui 

constraint . 

problen  is  to  ^axind^ 

_<_-t  us  con^  >e 

/  .nx  _  1,  .lv2    2/  ..  N\2 

g(o  )"p(d)+p  (8  )+...+( 

or 


g(en)  =  I  Pn(en)2  <F  , 

n=l 


where  p   is  constant  for  nth  stage.   1  tr< be  Ln 


follows 


where 


n  -  .n \ 2    /  n.r.i  ,    r.  .. 
P  v.b  )   =  (w  6  )(c  o  )  , 


w  =  weight  per  e- 

be  total  weight  of  ..- 

c  *  =  cost  per  elc. 

.e  cost  o:         e. 

Tr.is  gives  us  p   =  c  v 

To  apply  the  discrete  D  . 

follows ,  ] 

.n 

o   =  ituaber  o: 

re  dun- 


p  =  constant  of  nth  sta^e, 

x,  =  Amount  of  constraint  left  over  after  ailocati. 
first  n  stages, 

x_  =  Logarithm  of  the  reliability  upto  ar._  including  fia 
n  stages . 
This  N-stage  system  can  now  be  represented  by  the  following  performance 
equations. 


-  =  x£ A  -  pu(6")c,  n  -  1.  2,  ...  ,  N,  (U-22) 


n    ..-_    n/^.nN2 
x°  =  P, 

x2  =  x2_1  +  to  C1  -  t1-^)6*)  .  n  -  1,  2,  .  .  .  ,  H,   (U-23) 

x°  =  0. 


6  is  positive  integer,  i.e., 


Jn  >  1. 


The  objective  function  to  be  maximized  is  given  by 

S  =  I     An  (l  -  (l  -  Rn)6n)   =  4    . 
n=l 

This  gives  the  logarithm  of  maximum  reliability.   Com]  iring  s  juations 

(4-22)  and  (4-23)  with  the  equations  of  one  dimensional  proc given  by 

equations  (3-10)  and  (3-11),  we  obtain 

T(x^V)  =  x^1  -  Pn(en)2,  (.-2p) 

G(x^V)  -  to  (1  -  (1  -  R°     .  (14-26) 


— ■  -s 


ar(;    ••) 


3G( 


.. 


(x1  ,j    )   _,1 


1-  (1- 


[l*-i 


^^/  n~l  -rM 
3G(x,   ;o  ) 

— hn — "  °  • 

_  n-l 

3xi 


Let  us  define 


u  ■  1  -  BT  . 

.ng  equations  (^-27/  thrc  • 
of  one  dimensional  process  given 


-    '  ...    :-;  __ 

.. 
fi  -  -  ) 


-2p 


or 


(1  -    '    •  -    _    -  •••- 

CU*)6"  fa 


This  reduces  to 

„n+l    r-,        /l^+l\8i     •»  n  .     ..n+1       .        /,,nN6 

6  (1  -    (U        )  )  c_0 >  (1:.    u  \     (1  -   (IT )      1 

.    Qn+1  "    ^  n+lJ  ^        rTn  J     I  an     J    w      * 

(U^1)9  p      (uV 

Noting  that  all  terms  at  the  right  hand  side  are       if  9  is  known,  and 

we  need  to  find  6    ,  let 


n    „  TIn+l   .   /Tx-  3 
An  m    [JL,  rlr^,    j:-Ui  )  Qn 


» 


we  obtain 


Qn+1  (.         ,T.n+1.6n+\ 

o     [1  -  (U   )     I  _,  ,n 

(W1  "     ' 


en+1(i-  (un+1)e   )  =An(un+V 


or 


en+1  =  (An  +  6n+1)  (un+1)Qn+1  *  (h"- 

The  equation  is  solved  by  Newton's  method.   It  is  worth  not.     lat, 
while  8  ,  n  =  1,  2,  .  .  .  ,  N  are  in  reality  positive  integers  it  is  ass  —  d 
that  6  are  continuous  variables  in  obtaining  the  above  recur-- 
The  above  equation  may  be  written  as 

.n+1 


f(6n+1)  =  0n+1  -  (Aa  +  6n+1)   (Un+1)C"    ■  0  . 


As  LT    is  less  than  one,  it  can  be  easily  seen  that  f     )  v 

tonically  increasing  function  of  0*   .  =  0  w<  -s 

n+1 

negative.      As  we   increase  the  value  of  6        ,    this   function  will  assume   - 


.   klU6  , 

pass  t..  ro  one 

~xx   in  Fig  < 
app .     o  obtain  the  bo] 

Step  1.  Assume  6"1"  =  - , 
Step  2.  =1. 

Step  3.  Compute  k 

Step  k.  Compute 

f  {(o    ).}  ■  (8    ).  -  [A  t  [6 


ana 

(en~ 
f.{(rf»*)  >-l-  ((;r+1) 


Step  5«   Compute  a  new  trial  val.. 


2    -     _    / 


=  (6-     - 


Step  6.   Check  if 


"'2  "  <8    >ll  i 


(^  

ax 

is  sa 

Step  T.  The  value  of 

assumed  value  c     .  _,  us 

becomes  greater  than  '.- , 


^  .     .  .  .  ^         111 


•  I J .  ^ 


"J 


^> 


Fin    R 


wll  ^y.         .\.  >*    .  .    .  \^  •  .      »>  •  .  .  v*    .    .  .  \^  \_ 


for    solution    of     equation  [4-32), 


Step  6.       ~e 

■  p  -    I 

v  one  of  :.. 

is  (a)  go  to  step  9, 

stage,  go  to  b1 
Step  9.  Incr  i. 

Step  10.  Decrease  9  by  0.1 

again  greater  than     ,  when  a  go  to 

step  II, 
Step  II.  The  solution  has  reac. 

to  an  integer. 
A  computer  program  and  : 

ERICAL  EXAMPLES 

ree  stage  and  ei  .  s  are 

.od  developed  a:, 
oer  of  stages. 

...e  constants  assigned  for  a  three 

The  optimum  re 

e1  =  1.9000, 
e2  =  l, 

e3  =  z. 

Since  B  ,  n  ■  1,  2,  3      .  be  p. 
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Table  7.   Constants  assigned  for  3  stage  problem. 

^n 

n R g r__ 

1  .85  5.00 

2  .95  8.00  50 

3  .75  2.00 


•  =  2, 

»3-3. 

Thus  the  >ne 

from  each  of  the  above 

From  the  result  v 
used,  and  the  total  system  .  .    _ . 

j  constants  assigned  for  the  ej 
The  optimum  redundancy  obt 

e1  =  2.3000, 

62  =  1.57W, 

63  =  3.3109, 
6U  =  5.3310, 

e5  =  3.2^ 
e6  -  •  17, 
6T  =  3.15' 
6  =  2.53U1. 
Since  0  *,  n  =  1,  2,  .  .  .  , 

a1  -2, 
*2  =  z, 
e--  . 
•*-5. 


. 


Table   8.      Constants   assigned  for  8  stage  problem. 

n R^ ^ ?_ 

1  .85  5-00 

2  .95  8.00 

3  .75  2.00 
h                              .55                            1.00 

5  .60  7.00  300 

6  .65  3.00 

7  .70  U.00 

8  .80  5.00 


•'-3. 

>8-3. 

The  above  figures  being  - 
the  number  of  redunda; 
them. 

Here  we  note  that  267  Ui 
gives  the  optimum  relia: 


5.   OPTLMIZATION  OF  PROB]     IKVOLVI] 
SEPARABLE  FUNCTIC 


In  this  section  optimization  of  problems  with  a  separable  c     ive 

function  and  separable  and  monotonically  increasing  cons-:     .  is  con- 

sidered.  A  separable  and  monotonically  increasing  function  can  be  expressed 

as 

n 
f(x1,  x2,  .  .  .  ,  xQ)  =  I     t±(x.) 

i=l 

and 

f(x1,  x2,  .  .  .  ,  xi+l,  .  .  .  ,  xn)  >  f(xlS  x2,  .  .  .  ,  xQ) 

where  f . (x. )  is  only  a  function  of  x. .  Many  transformations  are  available 
to  reduce  the  function  to  a  separable  form.   In  previous  examples  ,  the 
function  of  total  system  reliability  is  reduced  to  a  separable  form  'oy   takin 
logarithm. 

A  general  problem  with  separable  functions  may  be  stated  as  folic  i  i 
Optimize 

s  =  I    fV),  (5-1) 

n=l 

subject  to 

N 

I    e?(en)  !*>•»    i  =  i,  2,  .  .  .  ,  s  .       (5-2) 

n=l  1       x 

There  are  various  methods  available  to  solve  such  types  of  problems,  -:.^..  as 
integer  programming,  the  gradient  method,  quadratic  pr  ip- 

proximation  etc.   In  this  section  the  discr  princ^  pplied 

to  such  problems . 


ve    i:.troauce    (s-  .  B ,    B    I 

........  .on.      ] 

D 

X.      = 

to  an..   til  -  •  •  , 

n 
x   ^,    -  tne   ace 

S+J. 

tht 
The  performance  equations   for  -st 

x^  =  x*"1  ♦  ^  a-  1,  2.   ....  *, 


I  -  1,  2, 


x°  =   0 


X.      <     D.  , 
1   —      i" 


x°       -  0. 

s+1 


The   objective   function  to  be  cptiz^izec 


s  =    V     .       - 

B»l 


Equation  (5-5)  gives 

c,    =  ,  2,    .    .    .    f  8, 

c  ^-    - 
The  Haailtonian  and  the   adjoi:. 


8+1 

T.n         r       n  n 

H     =     )      z.x. 

•    -,      i  i 


i=l 


=    I    z°  (x-1  ♦  el  (en))  ♦  <-+i  ««  +  ^8n»  ■         (5-7) 

i=l 


n  =  1,   2,    .    .    .    ,   N, 

n-1         311  n                      no                     tt                                       /■=:   fn 

z.       =  — —7  -  z. ,           n  =  1,  2,    .    .    .    ,  .. :                                   (?-a; 

1            «   n-1  1 
3x. 

1  i    -   1       2 

n-1  art  n  -,      o  »t  /  c   n    \ 

z„ ..   = =  z   ..  ,       n  ■  1,  2,    ,    ,    ,    ,  I.  (>-9a) 

s+1  n-1         s+1 

dX    , .. 
S+1 

From  equation   (5-6b) 

z*      =  1.  (5-9b) 

s+1 

Equations,  (5-9a)  and  (5-9b)  yield 

Zg+1  ■  1,  n  -  1,  2 N. 

This  reduces  equation  (5-7)  to 

Hn=  I    zMxJ"1  +  g?(en))  +x^  +  fn(en).        (5-10] 

1=1 

Differentiating  equation  (5-10)  with  respect  to  6  and  eq.u      to  ze] 
we  obtain 


n/  „n, 


8en     i=i  *   ae11 


or 


S    n  8 

1  ■■-   —  ♦ . 

I .e  Jth  co.. 

.e  end 
the  rest  of  then  are  . 

■  e±  ■      I  ■  1,  2,  ...  , 
J. 

From  equation  (5-o)  anc,  (;>-.    re  obtc. 

n  o 

i  *  J. 

..  «  1,  2 

Therefore,  equation  (5-11)  reduces  '.o 

n  3C,.(9  )         ) 

2   — ♦  ; ■  C.  -13] 

J 

The  procedure   for  solving  the  pre 

Step  1.     ..  og  some  v,. .  i  ,  we   obt  do 

fever,    equation    (5- 

Zj   =  Zjt  o  «  1,  2,   .    .    .    , 

Step  2.   Values  of  z     obtai.  -.,_,,.... 

from  equation  (5-- 
Step  3.   Check  to  see 

constraint  is  \ 

to  step  1.   If  the 

smaj.-i.er  value  for  - 


is  not  violated  and  some  other  cons        is  violat 
go  to  step  k.     If  the  active  constraint  reaches  its  1: 
none  of  the  other  constraints  are  violated  then  go  to  step  p. 
Step  k.      Replace  constraint  j  by  constraint  k,  i.e.,  now  con^":.-     -.  is 
active  and  the  rest  of  them  are  free.   'Therefore  in  equ_ 
and  steps  1,  2  and  3  replace  j  by  k  and  go  bach  to      1. 
Step  5.   We  have  reached  the  optimal  solution  with  the  constraint  in  cction 
as  the  only  active  constraint  on  the  syst 
In  the  above  iterative  process  we  nay  corr.e  across  a  case  where  con- 
sidering jth  constraint  in  action  violates  th^  constrain-  before  — 

jth  constraint  reaches  its  limit  and  considering  the  kth  constrain t    _..  .-ezic; 
violates  the  jth  constraint  before  the  kth  constraint  reaches  its  limit. 
This  indicates  that  the  optimum  lies  at  the  intersection  of  the  constraints 
j  and  k.   Under  such  condition  we  obtain  tie  following  equations. 
When  the  constraints  j  and  k  are  both  active  we  find  that 

zn  =  0,    i  =  1,  2,  .  .  .  ,  s  (5-lU) 

i  4   j,  k. 


This  reduces  equation  (5-11)  to 


n  3g,(Q  )    _  3gk(6  )    3f  (6  ) 

zn  —*- +  «?  -2L-_  + =  o  (5-15) 

j   30n       k    36* 
The  stepwise  procedure  to  solve  this  problem  La  given 

Step  1.   Assuming  some  value  for  6   in  equation  (5-15) «  •        -  •- 
relation  between  z  and 

Z]   =  hUk} 


1 


-    *«•  »    .    .    .    i 


and 


,  --,_>,... 


B   reduces   equation 


n   -  /  •• 

B  ■  li  2,  

i 


Suog;:l:  I  5-15  . 


)  . — ;  — 


=  0. 


2 
Step  2.  Assuming  son:e  value  for  6  in  eqv     I  (5-17 

n        2    n       ,   „ 
ana  z.     since  z.    =   z    ,   a  ■  1,  2.  •  .  .  ,  N. 

Step  3.  The  value  of  a.  obtained  la  

Q  , 

0  ,  n  =  3,  ■»,  .  .  •  , 
Step  k»      Check  if  constraint  J  or 

ditions  that  :.. 

(1)  r  of  the  two  CC  1, 

(2)  the  co:. 
ret. 

(3)  the  constrain 
re  ache. 


(k)     both  the  constraints  reac.  ■  _irr.it  s. 

If  condition  (l)  exists  then  30  to  step  5, 
if  condition  (2)  exists  then  go  "^°  step  £ 
if  condition  (3)  exists  then  go  to  step  7, 
if  condition  (U)  exists  then  go  to  step  8. 

Step  5.   Increment  9  to  some  higher  value  and  go  to  step  2. 
Step  6.   Increase  8  to  some  higher  value  and  go  to  step  1. 
Step  7.   Decrease  6  by  some  value  and  go  to  step  1. 
Step  8.   The  optimum  has  been  reached;  terminate  the  process. 

If  there  are  more  than  two  constraints  active  at  one  time,  the  above 
procedure  can  be  extended  to  deal  with  that  situation. 


6. 


..     . 


. 
constraints   are   alvayfl 
function   can  be  to 

system  reliability. 
objective   function 
(5-2)   and   (5-1)   respect! 
general   class   of  opti.  on  pre  I 

consider  a  prob.  two  stages  be 

illustrated  'zi~j  hand  computatioi 
..-jtage   system  and  can  be  ap~  . 
computer  program  is   given  in  Apr,.  -  . 

and  ten  stage   probler: 

^.tion  of  . 
redundant  units   increases    t] 
occupied.      L\       .    . 
straints  on  the  number  of  und 
presented  in  such  pre. 
Let 

C     =   Cost   of  one    .. 
of  on 

v"  =  Volume  o:' 

6     =  Number  o:'  e. 

n   constraints   can 


(1)  The  cost  constraint 

N 

g.  =  y  cnen  <  c  . 

X       n=l 

(2)  The  weight  constraint 

N 

v   n^n 

n=l 

(3)  The  volume  constraint 

g  =  f  vY<v. 

J   n=l 

The  problem  is  to  maximize  the  system  reliability  subject  to  the  above  con- 
straints. 

To  put  the  problem  into  the  discrete  maximum  principle  format,  let  us 
define  the  state  variables.   Let 

x,  =  Cost  of  elements  up  to  and  including  first  n  stages , 

x„  =  Amount  of  weight  allocated  to  first  n  stages, 

x„  =  Volume  occupied  at  the  end  of  first  n  stage, 

x,  =  Logarithm  of  reliability  up  to  and  including  first  n  sta 

The  problem  will  be  formulated  as  an  N  stage  problem.   The  performance 

equations  for  this  N  stage  problem  are  represented  as  follows 

n    n-1    n  n        ,   0         „  lc   .  , 

x=x    +cQ,    n=l,2,...,N,  --. 

xJ-O, 


..... 

x2  =  0, 

4  ■  °- 

x^  =  xj)"1  ♦  tn  (1  -  (i  -  B°        B-l.  2, 

x°-0. 

6  is  constrained  to  be  positive  intege 

o  ^  1. 

The  objective  function  to  be  :.. 

S  =  I     to  (1  -  (1  -  Rn)6  ) 

n=l 

Therefore  ve  obtain 

Bt  ■  0,    i  ■  1,  2, 

c.  =  - 
The  Haaiitoniaa  function  as 
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^  =  y  znxn 
.  ,   11 


i=l 


n    (  n-1   ,      n.n>i  n    r  n-1   ,      n„n^, 

=   z,     (x         +  c   8  j    +   z2    [x         +  v   o   J 


n  /-  n-1    nnn>    n  (  n-1      ,    ,    Dn,6\>  ^  _N 

+  z3  (x    +  v  d  J  +  z^  (x^   +  in   {1  -  (1  -  R  J   }J   ,       (o-T) 

n  =  1,  2,  .  .  .  ,  K, 


n-1    3H      n  , ,  Q , 

2i  =  71^=zi>  (o"°' 

3xi 


n-1    3H      n  ,,   _» 

Z2   =71^T=Z2' 
3x2 


n-1    3H      n  .,   , ^ . 

z3    "  7^1 =  v  (o-10) 

3x3 


n-1    Sn      n  ,  . 

zu  =  77^T=V  ^-^ 

3xu 

N 
zl+  =  ck   =  1*  (6-lla) 

From  equations  (6-11)  and  (6-lla)  we  obtain 

z£  =  1,     n  =  1,  2,  .  .  .  ,  ... 
Hence  equation  (6-7)  reduces  to 


1—2      n-\ 
*  I        ♦  i 

-12]  v. 
we  obtain 

=  0  =  z.c  t  ♦  -13 

1  -  (1  - 

Let  us  consider  now  bha 

rest  are  free.   Then  from  equation  [6-6] 


which  gives 


\  =   °2   ~~ 


*\   =  0,    n-  1,  2,  ...»  I, 


and 


.  c.  -  0, 


wniCw  t^ves 


z^  =  0,    a  -  1,  2,  ...  ,  .., 


refore  equation  (6-13        I  to 


-  .  -   :-         -  .: 

B  0  ■  T . 


Let  Un  =  1  -  Rn,   and  substituting  IT    in  eq  S-lU),  ve  obtain 


_  n  n  ,     r-(u    )        to  U  i  /  ^   -,  i     \ 

—  -  0  =  z,c  +    i j        .  (o-_-  a 

san  1  rs     fin 

39  1  -  (l/1)8 


Solving  equation  (6-lita)  for  z..  yields 


,      /TTns6      „   ITn 
n  1    flT  £r.L   «»  ,  e » 


c        1  -  (u") 


Solving  equation  (6-l4a)  for  6  ,  we  obtain 


n  n       .n 

zic      jn)e 


inU*  ±^n/ 


or 


n  n   «  „n 
z  c  +  SLnU 


(if)"  V 


or 


n  n 
z,  c 


ln( 1 } 

^  n  n   .  TTn  ' 


z  c  +  UnU 

en  =  — - 

ilnUn 


The  procedure  is  demonstrated  by  hand  computation     Lett —  for  a 


stage  problem.   The  constants  are  given  in  Table  9« 

Let  6  =1,  then  from  equation  (6-15),  we  obtain 


.  1  ' 

c     1- 


■  -    .  -    ■ 

9  l   1  -  (1  -  0. 


. 


From  equation  (6-c),  we  have 


1    2 
»1  "  Zl 


Therefore 


2  .    2. 

Zl  ~  "  c^ 


This  gives,  from  equation  (6-J 

in  (■ 


Z,  c 


2  2 

Z-C   +  • 


e2  = 


• 


«»t— 


=  1.57. 
Checking 
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Table  9.   Constants  assigned  for  2  stage  pro";.. 


Dn       n  n 

R       c       C      w  v 


1  0.9     9  5  1 

2  0.75    5      50     9      h$  1 


•  57 

=  .  . 


-  -       .      . 


>.03  -  1*5, 


-  -  -  1.57 


=  2.57  <  6. 


Thus  no  constraint  is  viola. 
61  =  2  . 


This  gives,  from  equation  (£-15), 


1   -2-3  _  2 
Zl  '    1    "1 


We  obtain,  froni  equation  (6-l6), 

62  =  3.3c  . 
Checking  if  any  constra. 


N 
x.  -9x2  +  5x3. 


=  3^.90  <  50, 


-  5  x  2  +  v  x   3. 


=  -....2  <  U5, 

=  2  +   . 


.o  constr-- 
This  gives 

:   -  ;• 
z-  = 


72 


Therefore 

e2  =  5.04  . 
Checking  if  any  constraint  is  violated 

3L  =  9  x  3  +  5  x  5.0lt 

=  1+2.2   <   50, 

x!J  =  5  x  3  +  9  x  5.0i* 

=  60.2  >  1+5, 
N 

x "  =  3  +  5.01+ 
3 

=  8.01+  >  6   . 

Thus  the  last  two  constraints  are  violated  but  not  the  one  cons .     I  to  he 

active.  So  we  make  this  constraint  free  and  assume  that  the  third  constraint 

is  active.  This  gives,  from  equation  {6-6) 

Zl  =  Cl  =  °> 

z"  =  0,    n  =  1,  2,  .  .  .  ,  N, 

N 
z2  =  c2  =  0, 

z^  =  0,    n  =  1,  2,  .  .  .  ,  N, 
and 

»^o  . 

Substituting  the  values  of  z  and  z0  into  equation 


z_  and  6  ,  we  obtain 


,   /,,nNG   .  Ttn 

v   l  -  (un)° 


and 


in   ( 


1        *        - 

-  ■ 


Nov  let  6=2 


3 


This  gives,  from  equation  (6-2 

62  =2.69  . 
Checking  if  any  constraint  is  viola;. 

x -  9  x  2  +  5  x  2.89 

=  32.45  <  50, 
x*«  5x2  +  9x2. fl 

-  36.1  <  45, 

=  2  +  2.89 

=  ^.o9  <  6  . 

6  be  3  which  gives 

2-   - 
Z3   "  999  ' 


We  obtain 


62  =  h.62    . 


7^ 


Substituting  in  the  constraints  ve  will  find  that  both  constraints 
two  are  violated  but  not  the  constraint  three.   Out  of  the  two  con     Qts 
violated  the  constraint  assumed  active  is  also  violated.  Therefore  to  know 
which  constraint  is  violated  first,  we  take  a  new  trial  value  for  8  less 
than  the  previous  value.  Let  6  =  2.5  which  gives 

1      2.3 
z~  =  - 


'3  "  "  311-5  ' 
and 

62  =  3.76  . 

Here  it  is  found  that  both  constraints  1  and  2  are  violated.  Let  the  new 

trial  value  of  9  be  2.4  which  gives 

1     2^3 
Z3  "  "  249   ' 

and 

e2  =  3-6  . 

None  of  the  constraints  are  violated  by  this  assumption  and  the  third  con- 
straint is  active.   Therefore,  the  optimum  policy  is 

61  =  2.4, 

e2  =  3-6  . 

It  should  be  noted  that  in  all  calculations  6  are  assumed  to  be  continu 
variables,  while  in  fact  they  are  positive  integers  only.   r_ ves 

61  =  2, 

e2  .  U. 

This  policy  is  not  violating  any  of  the  three  cc. 
solution. 


I 
However,  the  nu; 

systems  vith  an  arbitrary  Dumber 
constraints. 

The  constants  ass.  LO. 

The  optimum  redundancy  obtai..  La  L8  follows, 

81  =  2.9C 

82  -  •.  .20, 
e3  =  5.0> 

6U  =  3.7102, 

65  =  3.2U06. 
Since  6*  ,  n  =  1,  2,  .  .  .  ,5  should  be  posit 

i1 -3. 

e5  =  3. 

Tne  number  of  red;.  at   each 


From  the  result  we   ; 
..  9   units   of  siac/. , 
siaciK.   ana  volume   c. 
in   a  system  re_  1  of  0.9678. 
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Table  10.      Constants   assigned  for  5   stage  problem 


,n 


1 

.90 

8 

2 

.75 

9 

3 

.65 

o 

k 

.80 

7 

5 

.85 

8 

3  3 

1+  2 

155  9  132  -  50 

7  1 

7  2 


-  ■  . 
e2  =  - . 
e3  =  i.i»a 

o    =  2.6;:;.-, 

■  .  .. 

6      =   ^. 

6T   =   2.  31 

e8  =  U.   .     , 

-  •      32, 
u^u  =  3.3859. 

Since,    0,33  =  1,2,...,- 

6^3. 

*2  =  i. 

."-3. 

e5  =  2, 

e'  =  2, 

a6  =  5. 


Table  11.  Constants  assigned  for  10  stage  problem 

_n  n  _,        n  TT       n       „ 
n R c C w H v V__ 

1  .75  7  IT  1 

2  .60  h  19  1 

3  .95  U  25  1 
k  .80  10  Ik  1 

5  .85  k  175  15  ^75     1 

6  .90  8  16  1-0 

7  .80  6  21  1 

8  .55  3  11  1 

9  .70  6  18  1 
10  .65  2  lo  l 


e9  -  . 


B10 


The     tr  of  r  ay  be 

one  from  each  of  the  abc 

. 
.  of  1- 
vo^  ...its  with  a.  . 

reliability  of  O.ToTo. 


EXAMPLE  6.   MAXIMIZATION  OF  SYSTEM  RELIABILITY  . 
LINEAR  CONSTRAINTS 

In  this  example  the  method  described  in  section  5  •     -  applied  to 
the  problem  of  maximizing  system  reliability  subject  to  non- ar  con- 
straints.  The  method  developed  is  for  an  N-stage  syst^. -"..^s  are 

connected  in  series.  Each  stage  can  have  different  number  of  redundant 
units. 

In  the  last  example  we  assume  that  increase  in  cost  and  weight  due  to 
redundant  elements  is  directly  proportional  to  the  number  of  units  at  each 
stage,  that  is,  weight  and  cost  increases  linearly  with  the      r  of 
elements  at  each  stage.  This  is  not  necessarily  a  correct  assumption. 
the  number  of  units  is  increased  at  each  stage  it  requires  more  connect, 
equipment.  As  the  number  of  elements  at  each  stage  is  increased,  the  cost 
and  weight  may  increase  exponentially  due  to  the  additional  nook  up  required 
to  connect  these  elements.   Therefore  three  types  of  constraints  will  be 
considered. 
Let 

c  =  cost  per  element  at  nth  stage, 

w  =  weight  per  element  at  nth  stage, 

v  =  volume  per  element  at  nth  stage, 

8  =  number  of  elements  at  nth  staj 
(l)  The  constraint  of  the  type 

5     n     n  O 

g-,  =  I     P   (9  )  IP, 

x   n=l 

where  p  =  w  v  ,  the  product  of       per  u. 
nth  stage.   This  is  similar  to  one 


(2)  The  cost  con:; 


g?  =  I     cn(6n  *  exp. 

-..ere  c  0  will  be  the  cost  C 
the  additional  cojt  of 
(3)  The  weight  const 


g  =  I     wn6n  ex: 
3       n=l 

n  r 

where  w  6  *  will  be  the  « 

by  a  factor  exp(6  /U)  due  to  the  weight  Oj 

The  problem  is  to  naxinuz^ 
constraints.   To  formulate  I  .  discr 

we  define  state  variables  as  follov.  . 

Xl  =  Xl     ?      *    n  ■  1»  2,  .  .  .  , 

*[  l  P, 

n    n-1    n/  .n  ,   6  *,  /r  - 

X2  =  X2   +  c  ( 6  +  e   '  ) ,    r.  -_,_,..., 

0   r 
x2  =  0, 

'  <  c, 


x^  =  x^1  +  wVe^,    n  -  1,  2,  ...  ,  I,  (6-21, 

x3  =  °> 

xk  =  XS_1  +  Zn{1  "  (1  "  Rn)6*}  »    a  -  1,  2,  ...  ,  :;,    (6-22) 
x£  =  o. 
6  is  constrained  to  positive  integers  only,  that  is, 

en  >  l. 

The  objective  function  to  be  maximized  is 


S  =  I   An(l  -  (1  -  Rn)6  ) 
n=l 

k 

=  ^  =  2  ci^i  •  (6-23) 

i=l   x 

Therefore  we  obtain 

c.  =  0,    i  =  1,  2,  3,  .-24) 

ck   =  1. 
The  Hamiltonian  and  adjoint  variables  can  be  written  as 


r       n  n 

■  I  ■ 


♦  ■ 


♦  -.  ,-•-.       *  v  6  e    '  J  «*  +  ■ -     -   -  .j, 


3N  3 
n  -  1,  2,   3, 


z,         =  r  =    B.  ,  -26] 

oX. 


n-1        n 

zo     ■  — r~~  ■  zo> 

ox2 


,x3 


n-i 

oX. 


■     Cl|     =    lm 


Fro^  .  .    and   (    - 


-   1,  .  ,*.,...    , 


,rn  nf  n-i  J      n/Qn>.2>  rw  n-1         n,-n 

1^X1  P  J        Z2^X2       +  c    (6     +  e  )J 


nr  n-1   .      n,n  6  /h\  n-1  /•  ,  n 

+  z~[x,       +  w  6  e         J  +  Xi        +  fcn[l  -  (1  -  R   J      J    .  ^6-30, 

Differentiating  equation   (6-30)  with  respect  to   5*1  and  eqi  to   zero, 

we  obtain 


3h  _  n  n„n         n  n    (-        16   A> 

—  =  0  =  2z,p   9     +  z^c        1  +  —  e 
n  1  2        v  4  J 
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n  nr   6n/4  ,    1   0n   6n/^  -(1^)^(1-^  ,,   ..  , 

+   z_w     e  +  t-  6   e  +  .  (o-3U 

3      *■  4  '  „n 

,     / .    „n\ b 

1-U-R    ) 


Let     Un  =  1  -  Rn. 


Now  if  the   first  constraint  is   active   and  the   rest   of  them  are   free,  we    ...    in 
from  equation    (6-24)   the   following  relations. 


\  =  c±   =  0,  i  =  2,    3, 


which  gives 

\  -  0,  n  =  1,   2,    .    .    .    , 

i  =  2,   3. 

Therefore  equation    (6-31)   may  be  written  as 

nn 
/,Tn\D   «  ,,n 

i-(un)° 


Solving  equation  (6-32)  for  z.  yields 


1 

— 


2p   0        ,    , 

•  - 


Rearranging  the  ter 


V* 


b  A  ■  


Therefore  equation  (6-3h)  becc 


en  =  (uVV  t  An), 


or 


.. 


f(en)  =  en  -  .vr"/  -  An)  . 


The  equation  can  be  solved  as  exp_ 

If  the  second  cor., 
obtain  following  relations  fro;n  <_      .   - 


z\    =  c,  =  0,        L,  3, 


which  gives 


C  -  o,       -  l,  2,  ....  a 


i  -  1, 


Therefore  equation  (6-31)  may  be  written  as 


n  n 


z^c        1  +  t  e 


1      6n/k,  (-:nr   &r.b 


n,6r\   TTn 


l1  +  r e      J    = 


i-  u-v 


solving  equation  (6-3o)  for  zr   yields 


n 
Z2  = 


.n 


>(br-)5V" 


nf,    1  6  /^   .    /t.nj 
c  [1  +  ^  e    J   1  -  (b  ) 


Rearranging  the  terns  in  equation  (6-36),  we  obtain 


■36] 


(6-37) 


f(6n)  -  (1  +  i  e6^)  -  (bn/  ((1  ♦  i  e9^)  +  — -)   .    (6-3 

4  v     n  n  n  ' 


2 


The  equation  is  well  behaved  and  therefore  Newton's  method  car.  be  applied 

to  solve  it. 

Lastly  if  the  third  constraint  is  active  and  rest  of  then      ree 

we  obtain  from  equation  (6-2U)  the  following  relations 

N        «  ,   « 

z.  =  c.  =  0,     1=1,  2. 
l    l    '  ' 

This  gives 


z.  =  0, 

l 


n  =  1,  2, 
i  =  1,  2. 


Therefore  equation  (6-31)  may  be  written  as 


.. 


Solving  equatio:.        - 


■ ! ; —  .—  ,     ■ 

Rearranging  the  terns   in  eq I 

f(9")  =  ^   (i  ♦  i  en)  -  -  —   . 

3 

This  is  again  a  well  behaved  e 

The  computer  program  and  the  flow 

the  method  explained  in  section  5  are  gj  ...   

program  is  for  an  ii-sta^e  syst^ 

■ee  stage  and  five  stage 
for  the  three  stage  problem  are  give] 
obtained  is  as  follows  (see  «; 

e1  =  2.4000, 
e2  =  l. 

63  =  3.lk6$. 

Since  8',  n  =  1,  2,  3  should 


Table  12.      Constants   assigned  for  3  stage  probleci. 


Dn  n 

R  TO 


1 

.80 

1 

2 

.90 

3 

3 

.65 

k 

B8 


T  7 

TO  5  110  6  100 

9  6 


2 

=  2, 

The  number  of  redundant  elei 
one  from  eacr.  of  tr.e  above 

Froa  result  ve 
is  52  with  18  units  of  slac.-_. ,  the  cos*-  of  the  s; 
units  of  slacks  and  the  weight  of  th 
slacks.   This  policy  results  In  -  total  syi  . 

The  constants  assigned  for 
Tne  optimum  redundancy  obtained  is  as  t   Hows,  (see  Appe.. 

61  =  2.6000, 

62  =  2.28] 

u3  =  2.0075, 

e4  =  2. 

e5  =  3.3981. 

Since  6  ,  n  ■  lt  2,  .  .  .  ,5  she 

e2  =  2, 

oJ  =  2, 

.'-3. 

65  =  3. 


90 


Table  13.   Constants  assigned  for  5  stcj-  problem. 


Rn 


1 

.80 

1 

2 

.65 

2 

3 

.90 

3 

h 

.65 

h 

5 

.75 

2 

7  7 

7  5 

110      9  175      6      200 
9  6 

1*  9 


The  numbc.        .ndant 

. 
From  tht   result  ve 

and  volur..       -.-... 
.12  with  a  slack  of 

with  a  slack  of  7.52  US  . . 

of  0.90U5. 
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7.   PROPOSALS  FOR  FUTURE  STUDIES 

PROPOSAL  1. 

In  our  discussion  in  problems  1  to  6,  we  have  assumed  6"1  to  be  con- 
tinuous variables,  while  in  fact  obey  are  positive  integers  only.   Re  i 
off  to  the  nearest  integer  may  or  nay  not  give  the  true  optimum  value.   To 
avoid  this  ambiguity,  in  this  section  we  will  see  how  zero  one  integer  pro- 
gramming can  be  applied  for  rounc.  off.   Suppose  we  need  to  find  an  in 
solution  to  the  following  separable  programming  problem. 
Optimize 

f  =  I   fn(en),  (T-l) 

n=l 

subject  to 

I   S-(6n)  <  V    i  =  1,  2,  .  .  .  ,  r.      (7-2) 
n=l 

Let  the  integer  part  of  the  optimum  decision  obtained  by  solution  of  above 
problem  be  given  by  (6  )  ,   n  =  1,  2,  .  .  .  ,  N.   Let  9  denote  the  true 
optimum  integer  solution.   Then 


r  =  (en)1  +  (en)2  , 


where 


(en)2  =  0  or  1, 


The  new  problem  can  now  be  stated  as, 
Optimize 


F   = 

1 


subject,  to 


))  ♦  f  t     L  ♦  i)  - 


n=l  ^  a-1 


i  ■  1,  2 


ud 


(en)2  =  o  op 


Now  let 


n=l 


and 


.. 


P  and  b.  .  being  consta..      .e  above  for. 
Optiniize 


F2  «   I    ■     ' 
Ject  to 


I      (g"{(6n)1  +  1}  -  g^(6n)1}j(6n)2  <  b2i>    i  =  I,  2; 


n=l 


(7-6) 


and 


(en)2  =  0  or  1. 
Where 

F2  =  F  -  F1, 
and 

b„.  =  b.  -  b,.  . 

2i    i    li 

Both  the  functions  (7-5)  and  (7-6)  are  linear,  and  the  decision  variable  is 
restricted  to  zero  or  one.  Thus  this  problem  may  be  solved  by  zero  one 
integer  programming. 

PROPOSAL  2. 

In  the  application  of  the  discrete  maximum  principle  to  problems  wi«h 
separable  functions  as  given  in  section  5,  many  questions  need  to  be  answered, 
(i)  How  to  obtain  the  initial  estimate  for  6  , 
(ii)  If  more  than  one  constraint  is  violated  before  the  cons^raint 

considered  to  be  active  reaches  i*s  limit,  then  which  constraint 
should  be  selected  as  the  active  constraint  so  that  the  process 
of  convergence  to  optimum  is  fastest, 

(iii)  With  assumed  value  for  8  ,  if  no  constraint  is  vioia-ccc,  i is 

the  best  increase  in  assumed  value  of  6  ,  for  the  rapid  conver- 
gence to  the  optimum. 


PROPOw. 

possible  to  I 
parallel  redundancy  a] 

but  in  s*>. 

operation  when  the  or 

PROPOSAL  It. 

Optimal   choice   of  I    

that   at   any   stage   tfa  ne  poc. 

of  designs   that  are  possible  at  any  stage. 
designs   is  that  there  may  be  more    th  ype  of 

perforating  the  same  oper. 

and  volume,      Under  such   a  con. 
the  optimum  number  of  elemei 
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Appendix  A.      Con 

and  Resu- 


READ > 


PUT  :-:  -  o 

9  =C 

iUCViUi^ 


- 


n  =  2, . . . ,.;  froia 
equation   (^-10) 


COMPUTE  >;:•  - 

.;  -       I  w   (6  -- . 
"~~ 


,..-.     ........ 


"1? 

-     I  0.1 


^w-  ~-_ 


Table 


Program 

.  oi 


-ol 


R 

0 

X 

] 

Xl 

UR(I) 

-   - 

Rn 

IRED(I) 

6n 

-   1 

EEA(I) 

en 

M 

R(I) 

Rn 

w(i) 

w° 

ITHTA 

Rec 

of  e1. 

Wei^ 
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TABLE   A2.    COMPUTER  PROGRAM  FOR  E 


C   MAXIMIZATION  OF  SYS.  REL .  SU3  JECT  TO  /.  T  C 

DIMENSION  R( 10) tW(10) tX(10) iTHETAt 10) tUR(lO) »IRED(10) 

1  FORMAT(  15  ,F10.2) 

3  F0RMAT(F6.4,F5.2) 

104  FORMAT (28HDATA  FOR  THE  PROBLEM  FOLLOWS) 

105  F0RMAT(13HN0  OF  STAGES=I3) 

106  FORMAT (37HT0TAL  ALLOWABLE  WT  IN  REDUNDANT  ELE.=F6.2) 

100  FCRMAK1H  ) 

107  F0RMAT(8HoTAGE  N05X 1 1 HREL I  AS  I L I TY5X 1 1HWT/RED. ELE . ) 

108  F0RMAT(3XI3»9XF6.4,12XF;>.2  ) 

109  FORMAT (28H0PT I  MUM  POLICY  IS  AS  FOLLOWS) 

101  F0RMAT(6HSTAGE  »I3»3X21HNC  OF  REDUNC   IT  ELE.=  »F7.4»3X.3HUSE» I  5 ) 
103  FORMAT (9nT0TAL  WT= »2XF9.2 » 3X9HSYS.REL.= » 2XF6.4 ) 

REL=1. 

W  T  =  0 

M=0 

THETA( 1 1=0 

READ1 »N»PATT 

DO  2  I  =  1»N 

READ  3»Ri I ) »W( I ) 

2  CONTINUE 
PUNCH  104 
PUNCH  105,N 
PUNCH  106,PATT 
PUNCH  100 
PUNCH  107 

40  THETA( 1 )=THETA( 1 )+l. 

31  X (  1  )=PATT-W(  1 }*(THETA( 1 )-l. ) 

K  =  N-1 

UR(  1  )  =  1.-R( 1 ) 

DO  13  I=1»K 

UR(  1+1 )  =  1.-R{  1  +  1 ) 

Q=1.+IW(I)/W(I+1))*( LOG(UR( 1+1 ) )/LOG(UR( I ) ) ) 

P  =  1./UR(  I  )**THETA(  I )-l. 

THETA( I+1)=-(LCG(P*(Q-1.)+1. )/LCG(UR( 1+1 ) ) ) 

X(I+1)=X(  I  )-W(  I+1)*(THETA(  I  +  D-l.) 
13  CONTINUE 

I F ( X { N ) )  16, 19,15 

15  IF(M)  40,40,19 

16  THETAI 1 )=THETA( 1 )-.l 
M  =  M+1 

GO  TO  31 

19  DO  71  1=1 ,N 

71  PUNCH  10  8,  I  ,R(  I  )  ,W(  I  ) 

PUNCH  100 

PUNCH  109 

PUNCH  100 


52 

53 
54 
56 
51 


61 


61   1  =  1  i 
IT-iTA=T       [) 
DUMMY" ITH1 

. 
IF(DIFF-.5)  52.53.53 
! RFD( I )=DUMMY-1. 

1RED(  I  )  =. 

IF      M  I  )  )  56.56..  . 

rednt=thfta( i )-: . 

NCH  LOltltREDNTi      I  ) 

--IREu 

*EL*( 1«-(1«-R( I ) )**{RED+1«) ] 

P'oXCr 

punch   iC3.aT.rel 

END 


PLE  DATA  CARDS   ***»* 


N     PA~ 
5     104.0 


R 
.90 
.75 
.65 
.80 
.85 


W 

8. 

9. 

6. 

7. 

8. 


TABLE  A3. 


COMPUTER  OUTPUT  OF  FIVE  STAGE  PR.  . 


DATA  FOR  THE  PROBLEM  FOLLOWS 

NO  OF  STAGES=   5 

TOTAL  ALLOWABLE  WT  IN  REDUNDANT  ELE.=iC4.CC 


STAGE  NO 
1 
2 
3 

5 

RELIABILITY      WT/RED.ELE. 
•  9  0 C  C              8.00 
.7500              9.  CO 
.6500              6. CO 
.8000              7.00 
.8500              £.00 

OPTIMUM 

PC 

:LICY 

IS  AS  FOLLOWS 

STAGE 
STAGE 
STAGE 
STAGE 
STAGE 

1 
2 
3 
4 
5 

NO 
NO 
NO 
NO 
Nt: 

OF  REDUNDANT  FLE.= 
OF    REDUNDANT  ELE.= 
OF  REDUNDANT  ELE.= 
OF  REDUNDANT  ELE.= 
OF  REDUNDANT  ELE.= 

1  .60C 
2.8691 

4.230^ 
2.5o06 
2.0539 

USE 
USE 
USE 
USE 
US:l 

2 

- 
2 

TOTAL  WT= 

« 

LCA.Co    SYS.REL.= 

.9650 

TABLE  A4.    CCVPUT-      "PUT  CF  TEN  : 


Data  fcr   the   prc    .        FCllc 

nc  cf  stages=   1~ 
total  allc  .  . 


STAGE  N 

RE      .LITY 

1 

.  71 

2 

•  86 

- 

3 

. 

4 

. 

5 

. 

8. 

6 

. 

16. 

7 

. 

. 

8 

•  94  - 

:  l  . 

9 

■ 

. 

•  86 

cdt  : 

PC 

LICY 

FGLLCWI 

1 

CF    REDUND/ 

- 

. 

2 

. 

. 

27 

• 

ST* 

3 

' 

. 

= 

-  kGE 

NC 

- 

. 

SE 

5 

. 

=  2. 

2 

6 

. 

- 

3E 

7 

OF 

. 

- 

2 

STAGE 

0 

CF 

. 

= 

STAGE 

9 

. 

-  .  . 

STAGE 

IC 

CF  RE 

. 

= 

. 

TOTAL  ••"  = 

: 

198.      SYS.REL.= 

. 

3.06 


Appendix  B.   Computer  Flow  Chart,  Symbol  Table,  Computer  Pro 
and  Results  for  Example  2. 
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- 
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Table  Bl.  Symbol  table  for  computer  pr 


Program 
Symbol 


Mathematical 
Symbol 


-  vwiid. uior. 


REL 

R 
s 

WT 

PATT 

x°  =  W 

Ml 

m  -  1 

X(N) 

N 
Xl 

UR(I) 

1-Rn 

IRED(I) 

6n-  1 

THETA(I) 

en 

R(I) 

Rn 

w(i) 

n 
w 

ITHTA 

M 

Nl 


Total  system  reliability. 

Total  weight  of  the  system. 

Total  allowable  weight  in  redundant 

elements. 

Maximum  number  of  redundant  elements 
allowed  per  stage. 

Weight  left  over  after  allocating 
K  stages . 

Unreliability  of  one  element  at  nth 
stage. 

Redundant  elements  at  nth  stcge. 

Number  of  elements  at  nth  stage. 

Reliability  of  one  element  at  nth  stage. 

Weight  per  unit  at  nth  stage. 

Integer  portion  of  6*  . 

Dummy  variable  to  keep  track  in  variation 
of  a1. 

Number  of  stages  for  which  6n  is  less 
than  m. 


3LE  6. .    c:    :r  PR. 


C  N  Zr    SYS.REL.SUB. 

d;         .... 

1  FCRMATl  I5tF10.2) 

3  ....  2  ) 

:  3  ) 

1C5  FORMAT (  13m\C  OF  ST/ 

106  F  'CTAL  AL  6.2! 

100  FCRMA7     ) 

JoHSTAGL  N  .  > 

1C8  FORMAT ( 3X  :  ,     -  .  •  >       .2  ) 

109  FCRMAT(28HCF 

1C 1  FORMAT (6MSTAGE  » .   i  ... 

I  2  FCRMAT(9HTCTAL  WT= »2XF9.2 t 3X«  .   . 

L  =  l 
REL=1. 

r  a  ( l  > »  c 

.'  »P\TT 
RFAD4,'- 

r=c+.5 

do  2  i  =  :  .-. 

EAO  3»R< I ) »W( I ) 

2  CCNTI.V 
PL'.    .0^ 

punch  105#n 

punch  1c6»patt 

PUNCH  11 3, Ml 
PUNCH  1 

'  A  (  1  )  :  W  1 )  ♦  1 « 

31  X ( 1 )=PATT-W( 1 l»(THETA< 1 )-l. ) 
K*N  - 

l-LJ     33,33»32 

33  t-  r/wm+i. 

GO    TO    25 

32  URU )=1.-R( 1 ) 
DC    13    I=l,< 

t I+1)«1«-»R( 1*11 
0«(W( I )/W( 1  +  1) ) *(LCG  tUR(  )  ) 

P«     .         -.(I  )**TV,ETA(  I  )-l  . 

THFTA( 1+1 )*-(LCGCP*Q+l.)/LCG(UK  I      i) 

X (  1  +  1  ) =X(  I  )-W(  1  +  1 )*{ TH£TA(  i  +  1  )-l  . 
13    CONTINUE 
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N1=N-L+1 

IF(X(N1 ) )  16,19,15 

15  IF(M)  40,40,19 

16  THETA( 1 )=THETA( 1 )-.l 
MsM+i 

'   GO  TO  31 
19  11=1 

23  L1=N-L+1 

DO  21  1  =  11, LI 

IF(THETA( I )-C)  21*22,22 
22  IF( I-Ll )80, 35*35 
80  WET=W( I ) 

RE=R( I ) 

W(  I  )=W(L1 ) 

R(  I  )=R(L1 ) 

W(L1  )=WET 

R(L1  )=RE 

THETA( I )=THETA(L1 ) 
35  THETA(Ll) =C+.5 

PATT=PATT-W( Ll )*(THETA(Li)-l. ) 

L  =  L  +  1 

M  =  0 

IF(I-Ll)  10,21,21 
10  GO  TO  26 
21  CONTINUE 

1=0 
26  IF( I )  27,27,24 

24  11=1 
1=0 

GO  TO  23 
2  7  IF(M)  40,40,25 
2  5  DO  71  I=1,N 
71  PUNCH  108,1 ,R( I ) ,W( I ) 

PUNCH  100 

PUNCH  109 

PUNCH  100 

DO  61  1=1 . N 

ITHTA=THETA( I ) 

DUMMY=ITHTA 

DIFF=THETA( I >-DUMMY 

IF(DIFF-.5)  52,53,53 
52  IRED( I )=DUMMY-1. 

GO  TO  54 
5  3  IRED( I )=DUMMY 
54  IF(THETA( I ) )  56,56,51 
56  THETA( I )= 1 . 
51  THETA(I)=THETA(I)-1. 

PUNCH  101  ,1  ,THETA(  I  )  ,IRED(  I  ) 
RED=IRED( I ) 
WT=WT+W( I )*RED 
61  REL«REL*(1.-(1.-R(I) )#»(RED+1.) 

PUNCH  100 

PUNCH  103,WT,REL 

END 


^LF    D/  S       ***♦* 


75. 


R 

W 

. 

8. 

. 

9. 

.65 

7. 

ccnstraj  n   on 


***********  ************** 


TABLE    b3.  CCVPUTER    CUT 

WITh  RESTRICTED       - 

DATA    FOR    THE    PRO 
\~    OF    STAGES=       2 

AL^C'aABLE    WT  .=     7d. 

max  IDANT  rS   AuLCWEO   I 


STAGE  NC 

RELI     [TY 

•  S  E  . 

1 

. 

2 

. 

. 

3 

.6-: 

". 

CoT  ^CLICY     IS    AS    FCLU 


STAGE 

•3E 

STAGE 

1 
2 

TO' 

NC    OF    REOl  ... 

OF    REDU'.  -. 

72.  SYS.REl.= 
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table  b4.   computer  output  for  three  stage  problem 
with  theta  unconstrained 

data  for  the  problem  follows 

no  of  stages=  3 

total  allowable  wt  in  redundant  ele.=  75- 

maximum  redundant  elements  allowed  per  sta6e=  50 


STAGE  NO 

RELIABILITY 

WT/RED.ELE 

1 

.9500 

s.oo 

2 

.8000 

9.00 

3 

.6500 

7.00 

OPTIMUM  POLICY  IS  AS  FOLLOWS 


STAGE 
STAGE 
STAGE 

1 
2 

3 

TOTAL 

WT  = 

NO  OF  REDUNDANT  ELE.=  1.4000 

NO  OF  REDUNDALT  ELE.=  3.0085 

NO  OF  REDUNDANT  ELE.=  4.97o. 

70.0'j    SYS.REL.=  .9941 


USE 

1 

USE 

USE 

0 

TABlE  B5«    CC 

• TA  FOR  THE  PRC  .     FCLL: 
NC  OF  STAGES 

ALLOWABLE  ■  1  I 

4 
WT 

21. 

OPTIMUM  POLICY  IS  AS  FCLLC 

STAGE    1       CF  '-  .  =  ,. 


STAGE  NC 

RE. 

i 

•  7 

2 

.6; 

3 

. 

4 

. 

5 

. 

6 

. 

7 

.7 

8 

.9 

9 

.6 

10 

.7000 

STAGE 

2 

, 

.= 

2 

STAGE 

3 

CF 

- 

2 

STAGE 

4 

.  CF 

RE 

.  . 

■  -GE 

5 

.  ZF 

.  ..  n 

GE 

6 

ZF 

RE  . 

2 

•3E 

7 

ZF 

36 

3E 

c 

l\C  CF 

RE  . 

. 

2 

ST/ 

9 

ZF 

.' 

. 

-.. 

4 

-  -GE 

: 

ZF 

REDUI 

.:.= 

- 

TOTAL 

468, 

sys.rel.= 

.9' 
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TABLE  B6.    COMPUTER  OUTPUT  FOR  TEN  STAGE 
WITH  THETA  UNCONSTRAINED 


PROBLEM 


DATA  FOR  THE  PROBLEM  FOLLOWS 

NO  OF  STAGES=  10 

TOTAL  ALLOWABLE  WT  IN  REDUNDANT  ELE. =475.00 

MAXIMUM  REDUNDANT  ELEMENTS  ALLOWED  PER  STAGE= 


50 


STAGE  N 

RELIABILITY 

WT/RED.ELE. 

1 

.7600 

17.00 

2 

.8600 

19.00 

3 

.8800 

2  5.00 

4 

.7000 

14.00 

5 

.80^0 

15.00 

6 

.9000 

16.00 

7 

.  6  5  0  C 

21.00 

8 

.9400 

11.00 

9 

.79^0 

13.00 

10 

.8600 

18.00 

OPTIMUM 

PC 

:LICY 

IS  AS  FOLLOWS 

STAGE 

1 

NO 

OF  REDUNDANT 

ELE. 

=  3.CC 

USE 

3 

STAGE 

2 

NO 

OF  REDUNDANT 

ELE. 

=  2.1566 

USE 

2 

STAGE 

3 

NO 

OF  REDUNDANT 

ELE. 

=  1.6335 

5E 

2 

STAGE 

4 

NO 

OF  REDUNDANT 

EuE. 

=  4.0  014 

USE 

4 

STAGE 

5 

NO 

OF  REDUNDANT 

ELE. 

=  2.07o6 

USE 

3 

STAGE 

6 

NO 

OF  REDUNDANT 

ELE. 

=  1.8  3 

USE 

2 

STAGE 

7 

NO 

OF  REDUNDANT 

ELE. 

=  4.2207 

USE 

~r 

STAGE 

8 

NO 

OF  REDUNDANT 

ELE. 

=  1.5271 

USE 

2 

STAGE 

9 

NO 

OF  REDUNDANT 

ELE. 

=  2.8636 

USE 

3 

STAGE 

10 

NO 

OF  REDUNDANT 

FLE. 

=  2.1c 

USE 

2 

TOTAL  WT= 

468,00    SYS.REL.= 

.9782 

Appendix  C. 

and  Resu^ 


Mg.   CI.      Computer  1  low 


u - 


. ....  _  -  -  -   «* 


TGI.'.  cCI 


--  j.  j& » •  •  •  >« 

from  equation 

(-.21) 


iPUTE 
SYS.  


Table  CI.   .  ■  ol  table  . 


Program 
bol 


Dol 


Expl 


COST 

:) 

on 

.-. 

R(I) 

R° 

C(I) 

cn 

R 

s 

UR(I) 

1-  Rn 

ITHTA 

IRED(I) 

en-  i. 

. 


. 


io   be 

Reliability  c I 

sta 

j;er  pc.  - 
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TABLE  C2. 


COMPUTER  PROGRAM  FOR  EXAMPLE  3. 


C   C   MINIMIZATION  OF  THE  COST  OF  RED      IT  ELE       ..  .  _E  ACHIEVING 
C      A  MINIMUM  LEVEL  OF  RELIABILITY 

DIMENSION  R{ 10) tC(10) »THETA(10) »UR(10) » I  RED (10) 

1  FCRMAT(I3»F10.4) 

3  F0RMAT(F6.4*F5.2 ) 

103  F0RMAT(24HC0ST  OF  REDUNDANT  UN  I  TS  =  2XFlO  .  2  >  3X9HSYS  .  REL.  =2XF6  .  4  ) 

104  FORMAT (28HDATA  FOR  THE  PROBLEM  FOLLOWS ) 

105  F0RMATU3HN0  OF  STAGE5=I3) 

106  F0RMAT(26HMINIMUM  RELIABILITY  LEVEL=F6.4) 

100  FORMAT( 1H  ) 

107  FCRMAT(bHSTAGE  NC5X11HRELI AB I  LI TY5X13HCCST/RED.ELE. ) 

108  F0RMAT(3XI3,9XF6.4,12XF5.2) 

109  FGRMAT(28HGPTIMUM  POLICY  IS  AS  FOLLOWS) 

101  F0RMAT(6HSTAGE  »I3»3X21HNG  OF  REDUNDANT  ELE.  =  »F7.4»3X»3HUSE» 1 5 ) 
L  =  l 

COST=0 

M1  =  C 

M  =  0 

THETA( 1 )=0 

READ  l.NtRELB 

DO  2  1=1 »N 

READ  3,R( I ) ,C( I ) 

2  CONTINUE 
PUNCH  104 
PUNCH  105 »N 
PUNCH  1C6,RELB 
PUNCH  100 
PUNCH  107 

40  THETA( 1 )=THETA( 1 )+l. 

31  REL=1.-(1.-R(1) )**THETA<1) 
K=N-L 

UR(  1  )  =  1.-R( 1  ) 

32  DO  13  1=1, K 

UR(  1  +  1 )  =  1.-R(  1  +  1 ) 

Q=1.+(C(I)/C(I+1))*(L0G(UR(I+1))/L0G(UR(I))) 
P  =  1./UR(  I  )**THETA(  I )-l. 

THETA(  1  +  1  )=-(LCG(P*(Q-l. )  +  l.  ) /LOG( UR (  1+1 )  )  ) 
REL  =  REL*(  l.-(  l.-R(  1  +  1 )  )**THETA( 1  +  1  )  ) 
13  CONTINUE 

IF(M1 )  4,4,5 

4  IF(RELB-REL)  16,19,15 

15  IF(M)  40,40,19 

16  THETA( 1 )=THETA( 1 )-.l 
M  =  M+1 

GO  TO  31 
19  THETA(  1  )  =  THETA(  D  +  .l 
M 1  =  1 


5  DC  7: 

PUNCH  108.  I » R (  I  )  . C (  I  ) 

REL  =  . 

PUNCH  1 

PUNCH  I 

PUNCH  1 

0C  61  1  =  1  »i 

IThTA='       I) 

Dl     ITHTA 

DIFF«TH 

IFCDIFF-.5)  52»!  Ji 

52  IRET. .      ■y-i. 
GC  TC  54 

53  IREO( I JaOUMMY 

54  IF( THETA( I ) )  36»56»; . 
56  THETAt I )=1. 

51  THETA( I )=THETA( I )-: . 

PUNCH    lOl»ItTt  -  i ) 

RFD=IPrD( I ) 

CCST=CCST+C(  !  )*RED 
6 1  R  E I       ..-,,.-.:        + 1  •  ) ) 

PUNCH  100 

PUNCH  103.CCST.REL 

END 


SAMPLE  DATA  CARDS   ♦*»** 

N  LEVEL  OF  R 

R      C 
5     .90 

.65  2. 

.55 

.70  6. 

.75  7. 

.60  4. 
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TABLE  C3.    COMPUTER  OUTPUT  FOR  FIVE  STAGE  PROBLEM 

DATA  FOR  THE  PROBLEM  FOLLOWS 

NO  OF  STAGES=   5 

MINIMUM  RELIABILITY  LEVEL=  .9000 

STAGE  NO      RELIABILITY      CCST/RED.ELE. 

1  .6500  2.00 

2  .55-0  3.00 

3  .7000  6.00 

4  .75^0  7.00 

5  .6000  4.00 

OPTIMUM  POLICY  IS  AS  FOLLOWS 

STAGE    1    NO  OF  REDUNDANT  ELE.=  3.40C-    USE     3 

STAGE    2    NO  OF  REDUNDANT  ELE.=  3.9463    USE     4 

STAGE    3    NO  OF  REDUNDANT  ELE.=  2.0511    uoE     2 

STAGE    4    NO  OF  REDUNDANT  ELE.=  1.6406    USE     2 

STAGE    5    NO  OF  REDUNDANT  ELE.=  3.1501    USE     3 

COST  OF  REDUNLANT  UNITS=        56.00    SYS.REL.=    .902: 


tarle  C4.   :cmputer  •' 

data  for  the  pr: 

OF  STAGES=  ] 
MINIMUM  RELIABILITY  LEVEL=  . 


STAGE    N 

RE_ :      . _: ty 

CCST/"-     D.ELE. 

1 

. 

7.. 

2 

. 

)0 

3 

.    ■ 

<♦ 

. 

5 

.85 

<>• 

6 

. 

7 

. 

6. 

8 

•  55 

9 

.       00 

6. 

. 

•  65 

2.  - 

_ 

PC 

LICY 

IS    AS    FC 

STAGE 

: 

CF    Rl 

1. 

2 

GE 

2 

OF    Rl 

. 

GE 

3 

■ 

OF     R:      , 

EL. 

GE 

4 

Cr 

EL 

. 

OF    RE 

2 

6 

CF 

".= 

- 

STAGE 

7 

NC 

CF     REt 

2 

I     GE 

8 

• 

-. 

3E 

9 

JF     R: 

.       -7 

STAGE 

:  - 

CF 

L 

6 

CCST    ZF 

RE 

. :  ts= 

. 

. 

. 

. 
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Appendix  D.   Computer  Flow  Chart,  Symbol  Table,  Computer  Pre 
and  Results  for  Example  k. 


■ 


; 


- 


. 


•  - 


Table  Dl.   Symbol  table  for  computer  program. 


i-- 


Program 
Symbol 


Mathematical 
Symbol 


.  xplanation 


REL 

R 
s 

PR 

THETA(I) 

en 

PATT 

p 

R(I) 

Rn 

P(D 

n 
P 

X(I) 

n 
xl 

U(I) 

1  -  Rn 

A 

An 

THET 

f(en+1) 

THTA 

f (6n+1) 

ITHTA 

IRED(I) 

6n-l 

Total  system  reliability 

Total  constraint  limit  reach 

Number  of  elements  at  nth  stage 

Constraint  limit  allowed 

Reliability  per  element  at  nth  stage 

Constraint  constant  for  nth  stage. 

Constraint  limit  left  over  after  allocating 
n  stages 

Unreliability  per  element  at  nth  stage 

Constant  for  nth  stage  as  given  in 
equation  (^-32) 

Value  of  the  expression  given  by 
equation  (U-32) 

Value  of  the  expression  given  by  equation 
(U-32a) 

Integer  portion  of  6 

Number  of  redundant  elements  at  nth  staj  s 


TABLE  02.    COMPUTER  PROGRA 


C  riCN  CF  SYSTEM  REl 

»P(10)tX(10) i 

1  FORMAT ( 13tF6.2> 
3  FORMAT (•   .4, 

104  F:  ;ta  FOR  THE  PROBLEM  FOL. 

105  FORMAT( 13MN0  OF         ;  3) 

IAT( 18riC~  .2) 

:  •  ) 

107  FOR.'-        AGE  N05/  .  -  (  I  )  ) 

108  <VT(3XI3t9XF6.4tl2X  !  .. 

109  FC   \T(28H0PTIMUM  POLICY  IS  AS  FOL. 

I  :  FORMAT (6HSTAGE  .  I  3 » 3 X 2 1 r  .  .         ,15) 

103  F0RMAT(21HSLACK  IN  CONSTRAINT  ■■  . ■  .   L.«t2XF6.4] 

•  L  =  l. 
PR*C 

M 1  ■  ( 

ThETA( 1 >«0 

READ  1,N.PATT 
DO  2  1  =  1  ,N 

2  READ  3»R( I ) »P( I ) 
PUNCH  1 

PUNCH  105 »N 

PUNCH  1C6.PATT 

PUNCH  100 

PUNCH  1C7 
40  THETA(  1  )  =  THETA(  U+l, 
39  X (  1  )=PATT-P(  1 )*THETA( 1 )**2. 

<=N-1 

DO  13  1=1 .< 

U(  I  )=1.-R(  I ) 

UCI  +  1  )=1.-R(  1  +  1  ) 

A1=THETA( I )*?( I }*LOG(U( 1+1      L«-  .)) 

A2*P(I+1)*L0G(U( I ) )*U( I )*»T 

A  =  A  1  /  A  ? 
THE'     | )=1. 
37  T„rT  =  TnET       > - ( U (  I +  1  )**T-      I  +  1M*U 
THTA1«1«-U(I+1)»*THETA1  1+1 J 

THTA2*U( I+1)*«THETA( I+1)»LCG(U( !♦]  .  D+Al 

THTA=ThTAl-ThTA2 
D=THET/ThTA 
IF(ABS(D)-.C 1 )  36* 

35  THETAI  i  +  1  »«THETAI  I  +  D-l 
GO  TO  37 

36  X(I+1)*X(  I)-P(X+1)*THETA      **2. 
13  CONTINUE 

I  F  (  X  (  \  )  )  :  6  » 1 9  ♦  1  5 

15  IF(M1)  AC,.  ,19 

16  THETA( 1 )=TmETA( 1 )-. 
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Ml  =  l 

GO  TO  39 
19  DC  71  1  =  1  ♦  N 
71  PUNCH  106,1 ,R(  I  )  tP(  I  ) 

PUNCH  100 

PUNCH  109 

PUNCH  100 

DC  61  1  =  1  ,N 

ITHTA=THETA(  I  ) 

DUMMY=ITHTA 

DIFF=THETA( I ) -DUMMY 

IF(DIFF-.5)  52,53,53 

52  IRED( I )=DUMMY-1. 
GC  TC  54 

53  IRED( I ) =DUMMY 

54  IF(THETA(  I  ) )  56,56,51 
56  THETAI  I  )  =  1. 

51  PEDNT=THETA( I )-l. 

PUNCH  101,1 ,REDNT» IRED( I ) 

RED=IRED(  I  ) 

PR=PR+P( I)*(RED+1.)**2. 
61  REL=REL*(1.-(1.-R(I))**(RED+1.)) 

PR=PATT-PR 

PUNCH  100 

PUNCH  103,PR,REL 

END 


SAMPLE  DATA  CARDS  **** 


N   PATT 
3   50. 


R  P 

.85  5. 

.95  8. 

.75  2. 


TABLE  D3.    COMPUTER  OUTPUT  FOR  7-     ^TAGE  PR 


A  FCR  THE  PROBLEM  FCl 
NC  CF  STAGES=   3 
CONSTRAINT  LIMIT  - 

STAGE  NC  RELIABILITY  P(  I  ) 

1  .85 

2  .  9  5 

3  .75  2 . < 

odt:    POLICY  is  a?  follows 

5  OF  REDUN 

NO  OF  L  .  - 

:  OF  F  . -  . .        SE 

SlAC<  IN  CONST       =    ~.       SYS.f- .  !<►! 


stage 

1 

STAGE 

2 

STAGE 

3 

12£ 


TABLE    D4. 


COMPUTER    OUTPUT    FOR    EIGHT    STAGE    PROBLEM 


DATA    FOR    THE    PROBLEM    FOLl  OWS 
NO    OF    STAGES=       8 
CONSTRAINT    LIMIT    =300.00 


STAGE  NO 

RELIABILITY 

1 

.8500 

2 

.9500 

3 

.75^0 

4 

.55o0 

5 

.6000 

6 

.6500 

7 

.7000 

8 

.8000 

P(  I ) 

5.00 
8.00 

2.  CO 
1.00 
7.  CO 
3.00 
A. CO 
5.00 


OPTIMUM  POLICY  IS  AS  FOLLOWS 


STAGE 

i 

NO 

OF  REDUNDANT 

ELE.= 

1.300C 

USE 

1 

STAGE 

2 

NO 

OF  REDUNDANT 

ELE.= 

..748 

1 

STAGE 

3 

NO 

OF  REDU- 

ELE.= 

2.3109 

USE 

2 

STAGE 

4 

NO 

CE  REDUND- 

ELE.= 

4.3310 

4 

STAGE 

5 

NO 

OF  REDUNDANT 

ELE.= 

2.2511 

USE 

2 

STAGE 

6 

NO 

OF  REDUNDANT 

ELE.= 

2.63] 

USE 

3 

STAGE 

7 

NO 

OF  REDUNDANT 

ELE.= 

2.1578 

USE 

2 

STAGE 

8 

NO 

OF  REDUNDANT 

ELE.= 

1.5341 

USE 

2 

slack 

IN  0 

ONSTRAINT  =   13.0001 

SYS.REL.= 

.8384 

Appenoix  B«      Compute r 

and  Results  5« 


T34   .- 


tiy  )*i  .  - .  — 


_n     r. 

R    ,c    ,w   ,v   , 


;  ■•   ,   • 


compute  zr  fi  .. 


' 


D     — ^     —J. 


., 

COM  JTE 

. 

• 

_- 

'    .'  '-■- 

■ 

COMPUTE  :.". 

=     19' 


el-el-o.i 


. . .  ,:; 

i 


... 


Table  El.   Symbol 


Pro 

<ol 


.  cal 

t  ol 


-on 


R 
s 

:) 

en 

n 

?ATT(j) 

C(I,J) 

X(l,J) 

UR(I) 

1-   Rn 

L 

ITHTA 

IRED(I) 

en-  i 

,H1) 


Tot.- . 
.3er  of 
ber  of  coast] 

. 

— y  per  -~-ge 

Cor..  .  to  be 

i  Lon  of 

■ 
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TABLE   E2.   COMPUTER  PTCGRAM  FOR  EXAMPLE 


C   MAXIMIZATION  OF  SYSTLM  RELIABILITY  SUoJECT  TO  MULTIPLE  LINEAR  CONST. 
DIMENSION  THETA( 10) »PATT( 10 ) »R( 10) *C< 10*10 ) »K( 10 ) *UR( 10 ) *X( lOt 10) 

DIMENSION  IRED(IC) 
ICO  FORMAT ( 1H  ) 

101  FORMAT (41H0PT I  MUM  LIES  AT  THE  INTERSECTION  OF  CONST > I  3 ♦ 3HAND > I  3 ) 

102  FORMAT (28HDATA  FOR  THE  PROBLEM  FOLLOWS) 

103  F0RMATU3HNC    OF    STAGES=I3) 

104  F0RMAT(8HSTAGE  NO , 5X 1 1HREL I  AS  I L I T Y ) 

105  FCRMAT(3XI3»9XF6.4) 

106  FORMAT (28H0PT I  MUM  POLICY  IS  AS  FOLLOWS) 

107  FCRMATQ5X20HSYSTEM  RELIABILITY  =,2XF6.4) 

108  F0RMAT(6hSTAGE  »I3»3X21HN0  OF  REDUNDANT  ELE.  =  »F7.4»3X»3HUSE» 1 5  J 

1  FORMAT (2 15) 

3  F0RMAT(F8.2) 

5  F0RMAT(F6.4) 

7  F0RMAT(F10.2) 
L  =  l 

REL-1. 
M  =  0 

THETA( 1 )=0 
READ  1»N»N1 
DO  23I=i»Nl 
23  K( I )=G 

K(l  )=L 

DO  2  J=1,N1 

2  READ  3.PATTU) 
DO  6  I=1,N 
READ  5.R( I ) 

DO  6J=1.N1 

6  RFAD  7,C( I ,J) 
PUNCH  102 
PUNCH  103.N 
PUNCH  100 
PUNCH  104 

8  THETA( 1 )=THETA( 1 )+l. 
18  K2=0 

DO  9  J=1,N1 

9  X(1,J)=PATT(J)-C( 1*J)*(THETA( 1 ) ) 
UR(  1  )  =  1.-1M  1 ) 

Zl=UR(l)**THETA(l)*LuG(UR(l)) 
Z2=C(1 ,L) *( l.-UR(l )**THETA< 1 ) ) 
Z=Z1/Z2 

DO  11  I=2,N 

UR(  I  )  =  1.-R(  I  ) 

THT1=(Z*C(  I »L)  )/(Z*C(  I  •L)+LOG(UR(  I  )  )  ) 

THETA( I )=L0G(THT1 )/LOG(UR( I ) ) 

DO  11  J=1»N1 
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X{  i  ,J)=X(  I-l ,J J-C i  I  .  J)»l 

11  CONTINUE 

I F ( X ( N  *  L )  1    ... 

19  DC  12  I«li 

I F  ( X  ( N 1 1 )  1  1 3 1 1 2 1 

13  IF< I-<( I ) )  16.17,16 

16  < ( I  )  =  I 
L   . 

0 

THETA( 1)«THETA(1)-1, 

GO  TO  18 

17  k2=: 

12  CONTINUE 
IF(<2>  24,24,25 

2  4  I  F  (  X  )  6,3,. 

14  TnETA( 1 )=ThETA( 1)-. 1 
M  =  l 

GO  TO  18 

20  DO  71  1=1 ,N 

71  punch  1c5 #  i  »r( i  ) 
punch  ioo 

PUNCH  1 

PUNCH  1 

DO  61  1=1 ,N 

ITHTA»THETA( I ) 

DUMMY" ITHTi 

DIFF=THETA( I > -DUMMY 

IF(DIFF-.5>  52,53,53 

52  IRED( I )=DUMMY-1. 
GO  TO  54 

53  IRED( I »"OUMMY 

54  IF(THETA( I ) )  56,56,51 
56  THETA( I )=1. 

51  REDNT«THET«< I )-: . 

PUNCH  1C3,  I  ,REDNT, IRED(  I ) 
RED=IRED( 
61  REL=REL*(1.-(1.-R(I))**(RED+1.)) 
PUNCH  1C 
Punch  L07tREL 
GO  TO 

25  punch  ioi,k2»l 

30  STOP 
END 


SAMPLE  DATA  CARDS   ***** 


N 

5 


50. 

132. 

155. 


Nl 
3 


FIRST  CONSTRAINT  LIMIT 
SECOND  CONSTRAINT  LIMIT 
THIRD  CONSTRAINT  LIMIT 


DATA  STAGE  WISE 
.90 
3. 
5. 
8. 
.75 
2. 
4. 
9. 
.65 
4. 
9. 
6. 
.80 
1. 
7. 
7. 
.85 
2. 
7. 
8. 


TABLE  E3.    CC.VPuTER  OUTPUT  FOR  FIVE 


data  for  the  problem  follows 
nc  of  stages=  5 

stage  nc    reliability 

1  .9uw0 

2  .75 

3  .65 

5  .81, 

OPTIMUM  POLICY  IS  AS  FOLLOWS 

STAGE  1  NC  OF  RET  . -  L.                2 

STAGE  2  NO  OF  F  .611     USE 

STAGE  3  IC  OF  F  -    -          :E 

STAGE  4      .  OF  REP.  .  -  .  . 

STAGE  b  NC  OF  F  IT  ELE.=  2..-      USE     2 

SYSTEM  Rc'L.       •'  = 
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TABLE  E4. 


COMPUTER  OUTPUT  FOR  TEN  STAGE  PROBLEM 


DATA  FOR  THE  PROBLEM  FOLLOWS 
NO  OF  STAGES=  10 


STAGE  NO 

RELIABILITY 

1 

.7500 

2 

.6000 

3 

.9500 

4 

.8000 

5 

.8500 

6 

.9o„G 

7 

•  80o0 

8 

.55^0 

9 

.7000 

10 

.6500 

IPTIMUM  POLICY  IS  AS  FOLLOWS 


STAGE 

1 

NO 

OF 

REDUi.DANT 

ELE.= 

1. 

USE 

2 

STAGE 

2 

NO 

OF 

REDUNDANT 

ELE.= 

2.6784 

USE 

3 

STAGE 

3 

NO 

REDUNDA\T 

ELE.= 

.4220 

USE 

0 

STAGE 

4 

NO 

~F 

REDUNDANT 

EL£.= 

:  .6214 

US 

2 

STAGE 

5 

NO 

~F 

REDUNDANT 

ELE.= 

1.2735 

US 

: 

STAGE 

6 

NO 

~F 
\j  r 

REDUNDANT 

ELE.= 

.9286 

USE 

l 

STAGE 

7 

NO 

OF 

REDUNDANT 

ELE.= 

1.374 

US 

l 

STAGE 

8 

NO 

OF 

REDUNDANT 

ELE.= 

3.7186 

Li 

4 

STAGE 

9 

NO 

^F 

REDUNDANT 

ELE.= 

2.0632 

uSE 

2 

STAGE 

10 

NO 

~F 

REDUNDANT 

ELE.= 

2.3859 

USE 

2 

SYSTEM  RELIABILITY  = 


.7676 


Appendix  F.      Co::. 

and  Results 


F-     -      ti  n  r     ~ . 


r, 

„n     n     n     n 
,c    ,v    ,p 


-1 

*—  -L  i 


kn+1. 


1 


-.3.3. 


f  n+1 


'V+1] 


n+l 
BY   HEWTON'S 


j'Cl... 
1*1,2,3. 


..1 


,   ..-_ 


I 


Tab .    .   .   ol  tab 


Pro. 
Symbol 


RiX 

R 
s 

PR 

WT 

.'A(l) 

6n 

Nl 

R(I) 

RD 

U(I) 

1  -   R: 

P(I) 

n 
P 

C(I) 

n 
c 

w(i) 

n 

w 

X(I,J) 

a 

x, 

J 

z 

n 
z 

nus"! 

f(en) 

THTA 

r(en 

L 

;a 

IRED(I) 

1-   R1 

. 
Total  of 

e. 
Of  co:. 

/  per  unit  c    ... 
. 

Coefficic. 
Coefficic. 

over  c: 

-S. 

Ad.Join^  variable  a'. 

Value  v  - 

or   -    or  (U-  .  . 

;er  portion  o:* 


TABLE  F2 


COMPUTER  PROGRAM  FOR  EXAMPLE  6 


:   C   OPTIMIZATI 

DIMENSION 

DIMENSION 

ICO.  FORMAT  (  1H 

101  FORMAT(2I5 

102  F0RMAT(F1C 

103  FORMAT(28H 

104  F0RMAT(13H 

105  F0RMAT(8HS 

106  FORMATOXI 

107  F0RMAT(28H 

108  F0RMAT(6HS 

109  F0RMAT(15X 

110  FCRMATU1H 

111  FORMAT(9HT 
L  =  l 
REL=1. 

PR  =  0 
WT  =  0 
CST  =  0 
M  =  C 

THETA( 1 )=G 
READ  101, N 
DO  1000  1= 
1000  K( I )=0 
K(  1  )=L 
DO  1  1  =  1, N 

1  READ  102, R 
PUNCH  103 
PUNCH  104, 
PUNCH  100 
PUNCH  105 

2  THETA(1)=T 

3  K2  =  C 

X( 1,1 )=PAT 
X ( 1 ,2)=PAT 
X ( 1 ,3)=PAT 
U(l )=1.-R( 
DO  11  1=2, 
U( I )=1.-R( 
THETA( I )=1 
F=U( 1-1 )** 

4  GO  TO  (5,6 

5  Z=F/(2.*P( 
A=LOG(U(  I  ) 
THET=THETA 
THTA=1.-U( 


ON  OF  SYS  REL  SUBJECT  TO  NGN  LINEAR  MULTIPLE  CONSTRAINTS 

) »P(10) »W(10) fC(10) tX(lOtlO) tTHETA(lO) »U(10) 

( 10) ,K( 10) 


R(  10 

IRED 

) 

•  3F1 

.5,3 

DATA 

NO  0 

TAGE 

3,9X 

OPTI 

TAC.E 

20HS 

OPTI 

OTAL 


0.2) 
F5.2) 

FOR  T 
F  STAG 

NC5X 
F6.4) 
MUM  PO 

»I3»3 
YSTEM 
MUM  LI 

PR  =  F6 


HE  PROBLEM  FOLLOWS) 

ES=I3) 

11HRELIAEILITY) 


LICY  IS  AS  FOLLOWS) 

X21HNO  OF  REDUNDANT  ELE.=  , F7 . 4, 3X , 3HJSE , I  5 ) 

RELIABILITY  =,2XF6.4) 

ES  AT  THE  INTERSECTION  OF  CONST  * I3»3HAND» 13 

.2»9HTOTAL  WT=F6.2»10HTOTAL  CST=F6.2) 


»M»PATT1*PATT2»PATT3 
1.N1 


( I ) »P(  I  )  »C(  I  )  »W(  I  ) 

N 


HETA( 1 )+l. 

Tl-P( 1 )*THETA( 1)**2. 

T2-C( 1 )*(THETA( 1 )+EXP(THETA( 1 )/4.  )  ) 

T3-W( 1 )*THETA< 1 )*EXP(THETA< 1 )/4. ) 

1) 

N 

I  ) 

• 

THETA( 1-1 )*LOG(U( I  - 1 )  ) / ( 1 . -U ( I -1 )**THETA<  1-1 )  ) 

,7)  ,L 

1-1  )*THETA(  1-1  )  ) 

)/(2.*Z*P( I ) ) 

(I)-U( I )»*THETA( I )*(THETA(I)+A) 

I  )**THETA(  I  )-(THETA(  I  )+A)*U( I  )**TnETA(  I  )*LOG(U(  I  >  ) 


TO  8 


6    L*r/LC  *  .-!)*(  1.  -  ..  -  ,  )  ) 

_OG(U(  I  )  )/(2*C(  I  )  ) 
....  5*EXPl 
TH£T=B-U(  I  )**Tr 

THT1«(1«-     <  I  )**THETA( I )  )*EXP<7  .)/16. 

THT2=U( I )**1 
THTA  =  n-  72 

GO    TO    8 
7Z«F/(WCI-1)«EXP(TI  .  -       )  ) 

LCG(U( I))/(Z 
B«EXP(THETA( I)/4.)*(1«+.25*THETA( I) ) 
THET=B-U(  I  )*  - 
THT1**25*EXP(THETA( I )/4*)+B/4.- 

r2*U(i)«*THETA(I)*(«23*EXP(T  I       ■• 

TH1  '1-THT2 

8  D=THET/THTA 
IF(ABS(D)-.01)  iC.10,9 

9  THETA( I )=THETA( I )- 
GO  T~  • 

10  X(Itl)=X<  r-l»l)-P(I)*THETA( I )»*2. 

X (  I  ,  ?)  =X(  1-1,2 )-C(  I )*  ,  T  .  )  ) 

X(It3)=X( 1-1*3  -         7A(I)*EXP(       D/4.1 

11  CONTINUE 

I  F  (  X  (  N  ,  L  >  )  18,12.12 

12  DO  16  I=l,\l 

I F ( X ( N » I )  )  13,16,16 

13  I F (  I -K (  I  )  )  14,15.14 
1  A  K  <  I  )  «  ] 

L=I 

THETA(  1  )  =  THETA( 1  )-;  . 
GO  TC  3 

15  K2 

16  CONTINUE 
IF(K2)  17,17,27 

17  I F  <  M )  2,2,19 

18  THETA( i )=ThETA( 1 )-.l 
M=l 

GO  TO  3 

19  DO  20  1  =  1  ,. 

2C  PUNCH  106tl iR( I ) 

PUNCH  ICO 

PUNCH  107 

PUNCH  : 

DO  26  1=1 ,\ 

lTHTA  =  TnETA(  I  ) 

DUMMY* 1 ThTA 

DIFF«THETA( I ) -DUMMY 

IF(DIFF-.5)  21*22*22 
21  IRED( I »«DUMMY-1« 


1*2 


GO  TO  23 

22  IRED( I )=DUMMY 

23  IF(THETA( I ) )  24,24,25 

24  THETA( I )=1. 

25  REDNT=THETA( I )-l. 

PUNCH  108,1  ,REDNT»  IRED(  I  ) 
RED=IRED( I ) 

REL=REL*( l.-( l.-R( I ) )**(RED+1. ) ) 
PR=PR+P( I )*(REL+1.)**2. 
WT=WT+W( I)*(RED+l.)*EXP((RED+l.)/4.) 
2  6  CST=CST  +  C( I )*(  (RED+1. )+EXP(  (RED+1.  )/4. )  ) 
PUNCH  ICO 

PUNCH  109, REL 
PUNCH  111,PR,WT,CST 
GO  TO  28 

27  PUNCH  110,K2,L 

28  STOP 
END 


SAMPLE  DATA  CARDS   ***** 


N 

Nl 

PATT1 

PATT2 

PATT3 

5 

3 

110. 

175. 

200. 

R 

P 

C 

W 

.80 

1. 

7. 

7. 

.90 

3. 

5. 

8. 

.65 

4. 

9. 

6. 

.75 

2. 

4. 

9. 

.85 

2. 

7. 

P. 

TABLE  F3.    COMPUTER  OUTPUT  FOR  7.-.REE  S 


data  for  the  problem  foll: 
no  of  stages=  3 

stage  no    reliability 

1  .8000 

2 

3  .6500 

OPTIMUM  POLICY  IS  AS  FOLLOWS 

STAGE    1    -.0  OF  REDUNDANT  ElE.=  1.-         SE     1 

STAGE    2    NO  OF  r    _         .   .=     066 

STAGE         3         NO    OF    REDl  -.-    2.1465         US  2 

SYSTEM  REL ■       t    = 
TOTAL  PR=  52.00TOTAL  WT=  S7.57TCTAL  CST=  89.84 


1 


TABLE  F4.    COMPUTER  OUTPUT  FOR  FIVE  STAGE  PROBLEM 


data  for  the  problem  follows 
no  of  stages=  5 

stage  no    reliability 

1  .8GoO 

2  .8500 

3  .9C0C 

4  .7500 

5  .6500 

OPTIMUM  POLICY  IS  AS  FOLLOWS 

STAGE  1  NO  OF  REDUNDANT  ELE.= 

STAGE  2  NO  OF  REDUNDANT  ELF.= 

STAGE  3  NO  OF  REDUNDAMT  ELE.= 

STAGE  4  NO  OF  REDUNDANT  ELE.= 

STAGE  5  NO  OF  REDUNDANT  ELE.= 

SYSTEM  RELIABILITY  =    .9045 
TOTAL  PR=  83.00T0TAL  WT=192 .48T0TAL  CST=i46.12 


1.6C 

USE 

2 

1.2816 

USE 

1 

1.0075 

USE 

1 

1.6882 

USE 

2 

2.3981 

USE 

2 

OF. 


.    .  etrica] 

mbay,  ok 

3.E.    (Mechanical  EH 

- 
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ABSTRACT 

The  objective  of  this  report  is  to  obtain  solutions  to  reliability 
optimization  problems  subject  to  single  and  multiple,  linear  and  nonlinear 
constraints.  The  optimization  technique  employed  is  a  discrete  version  of 
the  maximum  principle. 

In  a  complex  system  where  a  number  of  stages  are  connected  in  series, 
each  stage  consisting  of  a  number  of  elements  in  parallel,  it  is  possible 
to  raise  the  reliability  of  the  system  to  one  by  employing  an  infinitely 
large  number  of  elements  at  each  stage.   However,  increasing  the  number  of 
elements  at  each  stage  increases  the  cost,  weight  and  volume  of  the  system. 
Since  resources  such  as  these  are  limited,  there  arises  a  systems  reli- 
ability optimization  problem.  The  two  types  of  this  problem  considered  in 
this  report  are  as  follows. 

(1)  Maximization  of  the  system  reliability  subject  to  (a)  a  s: 

linear  or  nonlinear  constraint,  (b)  multiple  linear  or  nonlinear  constraints. 

(2)  Minimization  of  cost  while  achieving  a  certain  minimum  level  of 
reliability. 

Examples  1,2,4,5  and  6  are  of  the  type  1.  In  Examples  1  and  2  a  si., 
linear  constraint  is  considered.  Example  2  is  an  extension  of  example  1, 
where  the  maximum  number  of  elements  that  can  be  used  at  any  one  stage  is 

restricted  to  some  value.  Example  h   considers  a  single  nonlinei        t. 

Examples  5  and  6  deal  with  multiple  linear  and  nonlinear  co..     nts 
respectively.   Example  3  is  of  type  2  where  a  linear  cost  function  is  min- 
imized while  achieving  a  minimum  level  of  reliability. 

In  all  these  problems  a  simple  computational  procedure  is  d 
based  on  the  discrete  version  of  the  maximum  principle.   Co: 


for  an  IBM  1620  ore  | 
developed  and  the  computer  pre 
Two  n'u. 


